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Abstract 

We study point and higher symmetries for the hydrodynamic-type systems with 
two independent variables t and x with and without explicit dependence of the equa- 
tions on t, x. We consider those systems which possess an infinite group of the hy- 
drodynamic symmetries, establish existence conditions for this property and, using 
it, derive linearizing transformations for these systems. The recursion operators for 
symmetries are obtained and used for constructing infinite series of exact solutions 
of the studied equations. Higher symmetries, i.e. the Lie-Backlund transformation 
groups, are also studied and the interrelation between the existence conditions for 
higher symmetries and recursion operators is established. More complete results are 
obtained for two-component systems, though n-component systems are also studied. 
In particular, we consider Hamiltonian and semi-Hamiltonian systems. 

1 Introduction 



In this paper we review our results of application of the group analysis of differential 
equations to systems of the hydrodynamic type. This concept includes a systematic 
study of point and higher symmetries, recursion operators, Hamiltonian structures and 
the usage of these results in order to obtain exact analytic solutions of these equations. 
Originally B.A. Dubrovin and S.P. Novikov meant by hydrodynamic-type systems quasi- 
linear systems of first order partial differential equations which possessed Hamiltonian 
structure fl],|2|]. Here we consider a more general class of hydrodynamic-type equations 
which specifically includes semi-Hamiltonian equations of S.P. Tsarev |3],f|] and explicitly 
dependent on t (or x) equations which are rich in symmetries. As a consequence of the 
last property they admit linearization and as soon as existence conditions for recursion 
operators are satisfied infinite series of their exact solutions can be obtained. Thus, they 
have as good integrability properties as the Hamiltonian equations. 

Hydrodynamic-type systems describe various physical phenomena: gas dynamics and 
hydrodynamics, magnetic hydrodynamics ||, nonlinear-elasticity and phase-transition 
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models ||, chromatography and electrophoresis equations from physical chemistry and 
biology |7|,|]. 

Applications of another kind are obtained by a representation of physically interesting 
higher order equations as integrability conditions of hydro dynamic- type systems: the Eu- 
ler and Poisson equations of nonlinear acoustics M, the Born-Infeld equation of nonlinear 
electrodynamics [flO[] , systems of relativistic-string equations Modern applications 

of hydrodynamic-type systems arise in the theory of averaging nonlinear soliton equa- 
lions i|]. 

An attractive feature of hydrodynamic-type systems is the possibility of a geometric 
formulation of their properties by close analogy with the Hamilton- Jacobi and eikonal 
equations of mechanics and optics. 

Symmetry group analysis of these systems leads naturally and algorithmically to asso- 
ciated differential-geometric structures: metric, connection, curvature, curvilinear orthog- 
onal coordinate systems and their transformations. If the Hamiltonian structure exists, 
then it turns out to be merely an aspect of this geometrical theory [|],[5],f|]. Thus, here we 
can see clearly a relation of transformation groups to geometry which S. Lie and F. Klein 
meant by their Erlangen program. 

A purely differential-geometric theory does not solve the problem of integrating 
hydrodynamic-type equations. We present here another approach to this problem based 
on a systematic study of higher symmetries and recursion operators |l2|-|l5l . The essence 
of this approach is the following. We pick out a class of those hydrodynamic-type systems 
which possess an infinite set of hydrodynamic symmetries depending on arbitrary solu- 
tions of a linear system of partial differential equations. Formulae for the corresponding 
invariant solutions determine a linearizing transformation, which reduces the problem of 
obtaining solutions of a nonlinear hydrodynamic-type system to a problem of solving a 
linear system. Then we may construct recursion operators which map symmetries again 
into symmetries if we assume in addition their existence conditions. They give rise to 
recursions of hydrodynamic symmetries and, as a consequence, lead to recursions for so- 
lutions of the corresponding linear system. Thus, we obtain the recursion formulae which 
allow us to multiply solutions of the original nonlinear hydrodynamic-type system. The 
action of recursion operators on the hydrodynamic symmetries explicitly dependent on 
x, t gives rise to infinite series of higher symmetries. It turns out to be possible to find 
the corresponding invariant solutions explicitly. Hence we obtain infinite series of exact 
solutions of the hydrodynamic-type system. They are analogous to the similarity solu- 
tions which are well-known in gas and fluid mechanics and therefore we may expect that 
they correspond to a physically interesting behavior. 

Symmetries, recursions, Hamiltonian structures and exact solutions of two-component 
spatially one- dimensional hydrodynamic-type systems were studied by the author mainly 
before S.P. Tsarev's publications [fl3] , [T6| p~^| . Symmetries, linearizing transformations and 
geometric theory of the multi-component hydrodynamic-type systems explicitly depen- 
dent on x, t have been constructed by the author in 1989 [ITU], EDI. Higher symmetries, 



higher conservation laws, their interrelations and the theory of recursion operators for 
multi-component hydrodynamic-type systems have been constructed by the author in 



1993, 1994 [pip,QapTJ. Finally in the recent article [22] we further developed the 



theory of integrability of the diagonal hydrodynamic-type systems with explicit t or x 
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dependence and presented a non-trivial example of such integrable system. There we also 
clarified and formulated precisely the concept of the hydrodynamic symmetries which are 
first order symmetries though being neither point, nor contact ones. 

In section 2 we study two-component diagonal hydrodynamic-type systems with ex- 
plicit dependence on t or x. We find an infinite group of their hydrodynamic symmetries 
and establish their existence conditions. By considering the corresponding invariant so- 
lutions we derive linearizing transformations for these systems. Second order symmetries 
and first and second order recursion operators are obtained together with their existence 
conditions which gives rise to infinite series of exact solutions of the studied systems. 

In section 3 the generalized gas-dynamics equations are studied which include in par- 
ticular the equations of the one-dimensional isoentropic gas dynamics. We obtain all their 
hydrodynamic and higher symmetries up to the third order inclusively. We construct first 
order recursion operators which lead naturally to the Lax representation of these equa- 
tions. The latter is used for obtaining the explicit formulae for infinite series of their 
invariant solutions. 

In section 4 a particular class of the two-component systems which may depend ex- 
plicitly on t and possess the Hamiltonian structure, the so-called separable Hamiltonian 
systems, is studied. We obtain their hydrodynamic and higher symmetries, linearizing 
transformation, second order recursion operator, the Lax representation and explicit for- 
mulas for the infinite series of exact (invariant) solutions. 

In section 5 we consider the n-component diagonal hydrodynamic-type systems with 
no explicit dependence on t and x which possess an infinite-dimensional group of the 
hydrodynamic symmetries, the so-called semi-Hamiltonian systems. We obtain their 
hydrodynamic symmetries together with their existence conditions and the associated 
differential-geometric structure. A linearizing transformation, first and second order re- 
cursion operators, higher symmetries and infinite series of exact solutions are constructed. 

Finally, in section 6 the n-component diagonal hydrodynamic systems with explicit 
dependence on t or x are studied. The hydrodynamic symmetries and linearizing trans- 
formations are obtained. 

2 Symmetries, recursions and invariant solutions 
for two-component hydrodynamic-type systems 

Here we study a two-component diagonal hydrodynamic-type system which may depend 
explicitly on time or space variable t or x. It is linear homogeneous in derivatives of the 
unknowns s(x,t),r(x,t) 



where the subscripts denote partial derivatives with respect to t and x. Functions and 
ip are real-valued and smooth and satisfy the following nondegeneracy conditions 





X 



(2.1) 



7^, 



<f) r (s,r,t)ip s (s,r,t) ^ 



(2.2) 
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Symmetries of system (|2 . 1| ) are generated by the Lie equations for one-parameter 
Lie-Backlund groups in a canonical representation p3, 24j 



tj — f ( T f q T V g(N) ~.(7V)\ 

1 t yy-^i u i °i 1 1 a n ' xi ■ ■ ■ i s x i 1 x )i 



(2.3) 



x T = t T = 0, (2.4) 

which are compatible with system ( [2.1|) . Here s x N ^ = d N s/dx N , = d N r/dx N and 

t is a parameter of a symmetry group: s = s(x, t,r),r = r(x, t, r). 

Compatibility conditions for systems (|2.1| ) and (|2.3|) s Tt = s tT ,r Tt = r tT take the form 
of the determining equations for characteristics (/, g) of symmetries 

D t [f] |(2.i) - <f>D x [f] - s x (<f>J + <f> r g) = 0, 
AMlQ - ipD x [g] - r x (ip s f + ^ r g) = 0, 

where D t and D x are operators of the total derivatives with respect to t and x 



(2.5) 



D = - + s - + r ^ + y(s^^- + r^^-] 

X dx +Sx ds + rx dr + {** ds^ + x dr^ ) ' 

A|(2.i)= |r + <k|- +i/jr x ^-+ (2.6) 
at os or 

+ £ ( ^[(ps^—TTT + D^[ifjr x 



and _D t is calculated with the use of system ( |2.1|) . 

2.1 Hydrodynamic symmetries of diagonal systems with 
explicit temporal dependence 

For hydrodynamic symmetries put N = 1 in the Lie equations ( |2.3| ) 

S T f(%: S, T, Sa;, ^rc); )V g(%: '"j 'Sx; ?"x); (2.7) 

and in determining equations (|2.5|) . Analysis of equations ( |2.5|) produces the following 
results. Define the functions 

$ r (s,r,t) = <f> r (s,r,t)/(<f>-i>), e,(s,r,t)=Vr(s,r,t)/(V-0), (2.8) 

(2.9) 



$(s, r, i) = 6(s)$ a (s, r, t) + d(r)$ r (s, r, i) + $ (s) 
9(s, r, t) = b(s)S s (s, r, t) + d(r)9 r (s, r, £) + o (r) 



(2.10) 



0(s, r, t) = 6(s)0 5 (s, r, i) + d(r)(f) r (s, r, t), 
$(s, r, t) = b(s)ip s (s, r, t) + d(r)ij) r (s, r, t). 

Definition 2.1. We call system ( |2.1|) a generic system (with respect to hydrodynamic 
symmetries) if its coefficients 0, ip do not satisfy the constraints 



t = /3(t)0 2 + e(t)0 + A(t), V* = + £ W + X(t) (2.11) 
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with arbitrary smooth functions /3(t), e(t), X(t). 

Results of the group analysis are formulated as the following basic theorem [TJ 
Theorem 2.1. A diagonal two-component generic hydrodynamic-type system ( |2 . 1| ) 
which may explicitly depend on t possesses an infinite set of hydrodynamic symmetries 
with a functional arbitrariness iff the following two conditions are satisfied: 
1. Coefficients <f>,i() of system (|2.1|) satisfy the equalities 



rt 



e 



St 



Ms 



[2.12) 



where (3 is arbitrary (real) constant, the functions $(s, r, t), 9(s, r, t) are defined by equa- 
tions ( |2.8| ) and partial derivatives with respect to t are taken at constant values of s and 

r; 

2. There exist four functions of one variable b(s), d(r), <&o(s), &o(r) which satisfy the 
equalities 

$ r = $ r ($ - e), e s = e s (e-<l) (2.i3) 

where the functions $(s, r, t), B(s, r, t) are defined by the formulae 
These symmetries are generated by the Lie equations 



s T = f = <j)(x, t, s, r)s x + b(s), 
r T = g = $(x, t, s, r)r x + d(r), 

where functions 0, ip are defined by the following formulae 
if / 



(2.14) 



4>(x, t, s, r) = a(s, r) exp j/? 
ijj(x, t, s, r) = c(s, r) exp j/? 
and if (5 = 



x + J 0(s, r, t)c?t 
o 

x + J ^(s, r, t)c?t 
o 



$(s,r,t), 



(2.15) 



+ -e(s,r,t), 



(x,t,s,r)= a(s,r) + J <ft(s,r,t)dt— 



x + J <fi(s, r, t)dt 
o 



ip(x,t,s,r) 



-$(s,r 
c(s, r) + J ip(s,r,t)dt— 



(2.16) 



-9(s,r) 



x + / ^(s, r, t)(it 
o 



Here the integrals with respect to t are taken at constant values of s and r and the 
functions 0,-0 are defined by the formulae (|2.10|) . The functions a(s, r), c(s, r) constitute 
arbitrary smooth solution of the linear system of equations 



a r (s,r) — $ r (s, r, 0) (a — c) , 
c s (s,r) = 6 s (s,r, 0)(c - a). 



(2.17) 



Remark 2.1. We can use the freedom in the definition (|2.8|) of the functions $, to 
transform equations ( |2.12| ) in the condition 1 of Theorem 2.1 to a more simple form 



$t(s, r, t) = p(p(s, r, t), Q t (s, r, t) = /30(s, r, t). 



(2.18) 
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Remark 2.2. A solution manifold of linear system (|2.17|) is locally parametrized by two 
arbitrary functions c\(s), c^ir) of one variable. They determine a functional arbitrariness 
in the definitions ( |2.15| ), (|2.16|) of hydrodynamic symmetries ( |2.14| ) for system ( |2 . 1| ) . 

System (|2.17|) always has a trivial solution a(s, r) = c(s, r) = c = const. 

The condition 2 of Theorem 2.1 always has a trivial solution 



b( s ) = d{r) = 0, 

$ = e 



$o(s) = e„(r) 

C , = $ -- 



Co = const, 



0. 



(2.19) 



Corollary 2.1. The condition 1 of Theorem 2.1 is necessary and sufficient for system 
J]) to possess an infinite set of hydrodynamic symmetries, generated by the Lie equations 



(x,t,s,r)s a 



r T = i[)(x, t, s, r)r a 



(2.20) 



which are linear homogeneous in derivatives. The coefficients 0, ip of these equations are 
determined by the following formulae 
if/3^0 



cj)(x, t, s, r) = a(s, r) exp < (3 
ip(x, t, s, r) = c(s, r) exp < (3 



x + J 0(s, r, t)dt 



x + J i[)(s, r, t)dt 
o 



+ c , 
+ c , 



(2.21) 



and if (5 = 



{x,t,s,r) 



a(s,r) + c 
4>(x,t, s,r) = c(s,r) + c 



x + j 4>(s, r, t)dt 
o 
t 

x + j i[)(s, r, t)dt 
o 



(2.22) 



Remark 2.3. The condition 1 of Theorem 2.1 with f3 = is met in particular for 
system 2.1 with the coefficients 0(s, r), ip(s, r) explicitly independent of t. Such a system 
always has an infinite set of hydrodynamic symmetries with a functional arbitrariness. In 
this case coefficients 2.22 of Lie equations have the form 



<fi(x, t, s, r) = a(s, r) + c [x + t(j)(s, r)) , 
ip(x, t, s, r) = c(s, r) + c [x + tip(s, r)] . 



(2.23) 



Remark 2.4- The condition 2 of Theorem 2.1 is that additional constraint which 
provides the existence of symmetries with the Lie equations (|2.14j ) linear inhomogeneous 
in derivatives. Every nontrivial solution of equations (|2.13|) generates such symmetries. 



2.2 Infinite-dimensional Lie algebra of hydrodynamic symme- 
tries and recursions of symmetries 



Let the system fl2.1|) possess two one-parameter symmetry groups generated by Lie equa- 
tions ( |2.14| ) with a parameter r and by Lie equations of the same form with a parameter 

f 

Sf = f = <j)(x,t,s,r)s x + b_(s), , 2 24 ^ 

r f = g = ijj(x,t, s,r)r x + d(r). 
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Here the coefficients <fr, ip are determined by analogy with the formulae (|2.15 ) or ( |2.16| ) 
if / 



(p(x, t, s, r) = a(s, r) exp |/5 
ip(x, t, s, r) = c(s, r) exp j/3 
and if (5 = 



t 

a; + / <p(s, r, t)dt 
o 
t 

x + J ip(s, r, t)dt 
o 



+ ~e(s,r,t), 



(2.25) 



(x,t,s,r)= a(s,r) + f <p(s,r,t)dt— 



x + J 0(s, r, t)dt 
o 



-$(s,r 

ijj(x,t, s,r) = c(s,r) + f ip(s,r,t)dt— 
— Q(s, r 

The functions S>, G are defined by analogy with the formulae ( |2.9| 



(2.26) 



x + J ip(s, r, t)dt 
o 



$ = 6(s)$ s + rf(r)$ r + $ (s), 
6 = 6(s)0 s + J(r)e r + o (r), 

and the functions 0, ip are defined by analogy with the formulae (|2.10|) 

(f> = b(s)(f) s + d(r)(f) r , $ = b(s)ip s + d(r)tp r . 
The functions G must satisfy the condition 2 of the theorem 2.1 

<S r = $ r (i - <§>), e s = e s (e-i>) 



(2.27) 



(2.28) 



(2.29) 



The functions a(s , r) , c(s , r) constitute arbitrary smooth solution of the same linear 
system ( |2.17| ) 

a r (s,r) = $ r (s,r,0)(a- c), , 2 3Q . 

c s (s,r) = 6 s (s,r, 0)(c- a). 

Let a = (/,<?) and a = (/, <?) be characteristics of the symmetries ( |2. 14j ) and ( |2.24| ). 
A canonical Lie-Backlund symmetry operator with the characteristic a is defined as fol- 



lows 24 



d , _ 9r „ nN d ,,,,, „ 9 



••• + W[/])i + (^D- " 



<9r. 



+ . . . 



(2.31) 



+ . . . 



Here operator of the total derivative with respect to t is calculated with the use of 
equations Q2.ll) . The formula for is obtained by a substitution of (/, g) for (/, g) to 
the formula ( |2.31| ). 
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The usual Lie commutator 



[Xgj X a 



Xt, 



X, 



(2.32) 



generates the commutator of symmetry characteristics 



a = {a, <r\ = a [a\ - a [erj. 
Here a' denotes the operator of the Frechet derivative |]25|| ) 

00 da n 



^'U = Y,^m D i (a, 13 = 1, 2), 



3=0 



(2.33) 



(2.34) 



where <j\ = /, 02 = g,u l = s,u 2 = r. The commutator ( |2.33| ) of symmetry characteristics 
for system (|2.1| ) is again a characteristic a = (f,g) of some symmetry for this system 
generated by the Lie equations with the parameter f: Sf = f,Tf = g where 



/=[ct,(tJi= /- J-T+9-, 9~r + 

as os or or 

- df df 

g=\a,a\ 2 = f-= f-= V g-= g-= h 

os os or or 

+D -®wr D '®w> 



[2.35) 



(2.36) 



Theorem 2.2. (see [TJ). 

A commutator of hydrodynamic symmetries (2.24) and ( |2.14 ) with the characteristics 
a and a is again a hydrodynamic symmetry of system ( |2.1fl with the characteristic a = 
[a, a] = (f,g), which is generated by Lie equations 



Sf = f = 4>(x, t, s, r)s x + b(s), 
Tf = g = if)(x, t, s, r)r x + d(r) 

where the coefficients are determined by the following formulae 
if (5 ± 



<f)(x,t, s,r) = a(s,r) exp j/3 
+-4(s,r,t), 

i)(x,t, s,r) = c(s, r) exp j/3 
+-Q(s,r,t), 



x + J <fi(s, r, t)dt 




x + J ^(s, r, 




(2.37) 



(2.38) 
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and if (3 — 



t % 



(x,t,s,r) = a(s, r) + J r, t)dt— 



-$(s,r) 



x + J (j)(s, r, t)<it 



(2.39) 



il>(x,t, s, r) = c(s,r) + J ip(s,r,t)dt— 



-Q(s,r) 



x + J ip(s, r, t)dt 
o 



and for any value of /3 

5(*) = 

Here the following notation is used 



b(s) b(s) 
b'{s) b'{s) 



d{r) 



d(r) d(r) 
d'{r) d'(r) 



(2.40) 



a s, r 



c s, r 



0(s, r, t) = fo(s)0 s (s, r, t) + d{r)(f) r (s, r, t), 
tjj(s, r, t) = b(s)i[) s (s, r, t) + d(r)i[) r (s, r, t). 

<l(s, r, t) = b(s)& s + (J(r)$ r - b(s)& s - rf(r)l> r , 
9(s, r, t) = b(s)<d s + d(r)O r - b(s)Q s - d(r)Q r , 

b(s)[a s (s,r) - $ 8 a(s,r)] - b(s)[a s (s,r) - ® s a(s,r)}+ 
+$ r [d(r)c(s, r) - d(r)c(s, r)] + $ (s)a(s, r) - $ (s)a(s, r) 



(2.41) 



(2.42) 



(2.43) 



^( r )[ c r(s, r) — r c(s, r)] — d(r)[c r (s,r) — r c(s,r)] + 
+6 s [6(s)a(s, r) - b(s)a(s, r)] + o (r)c(s, r) - o (r)c(s, r), 

$ = $( Sj r , 0), = 0(s, r, 0) if ^ and $ = $(s, r), = 0(s, r) if /? = 0. 

Corollary 2.2. The formulae Q2.42|) and Q2.43j ) determine recursions of solutions for 
the linear systems fl2.13|) and (|2.17|) respectively. They map any pair of solutions of the 
corresponding system again into its solution. 

Consider now a special case 



b( s ) = d{r) = 0, $ = = c = 0, (b{s), d(r)) ^ (0, 0). 



(2.44) 



Corollary 2.3. Let the condition 1 of the theorem 2.1 be satisfied for the system 
( |2.1D and inhomogeneous in derivatives hydrodynamic symmetry of the form ( |2.24j ) exist. 
Then the Lie commutator of the symmetry Q2.241 ) and a homogeneous in derivatives 
hydrodynamic symmetry of the form ( 2.20|) -( 2.22|) with Co = is again a homogeneous 
symmetry of the same form 



Sf = f = (f>(x,t, s,r)s, 
r f = g = tp(x,t,s,r)r 



(2.45) 
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where for any real value of (3 



4>{x, t, s, r) = a(s, r) exp < (3 
V>(x, t, s, r) = c(s, r) exp < (3 



x + j 4>(s, r, t)dt 



x + J if)(s, r, t)dt 
o 



;2.46) 



and 



a(s, r) 



c{s, r 



b{s)a s (s, r) - [6(s)$ a (s, r, 0) + $ (s)]a(s, r)- 
— d(r)<& r (s, r, 0)c(s, r), 

d^ev^r) - [d(r)6 r (s,r, 0) + 9 (r)]c(s, r)— 
-6(s)e a (s,r,0)a(s,r), 



;2.47) 



Thus, the Lie commutator with the inhomogeneous in derivatives symmetry (|2.24j) 
is a linear operator acting on the space of homogeneous hydrodynamic symmetries with 
Co = 0. It generates the recursion ( |2.47|) of solutions for the linear system ( 2.12|) . 

Remark 2. 5. Inhomogeneous in derivatives hydrodynamic symmetry generates via Lie 
commutator a recursion operator for homogeneous hydrodynamic symmetries. 

Corollary 2.4. If the system ( p.l|) satisfies the condition 1 of the theorem 2.1, 
then it possesses an infinite set of mutually commuting hydrodynamic symmetries with a 
functional arbitrariness which are homogeneous in derivatives with cq = and have the 
form ( ggg ), ( ggp - 

In particular, if for (3 = the coefficients (f)(s,r),ip(s,r) of the system (|2.1| ) do not 
depend explicitly upon t, then we reproduce Tsarev's result about the commutability 
of hydrodynamic-type flows without explicit dependence on t, x. 



2.3 Hydrodynamic symmetries of diagonal systems with 
explicit spatial dependence 

Consider a diagonal two-component hydrodynamic-type system with explicit dependence 
on x 

s t = r, x)s x , r t = i[>*(s, r, x)r x (2.48) 
Define functions $(s, r, x), 0(s, r, x) by the following equations 

6 s (s, r, x) = r, x)] s /[^* _1 (s, r, x) - r, x)] 

Theorem 2.3. A diagonal generic hydrodynamic-type system (|2.48[) which may ex- 
plicitly depend on a space coordinate x possesses an infinite set of hydrodynamic sym- 
metries with a functional arbitrariness locally parametrized by two arbitrary functions 
Ci(s), c>z{r) of one variable iff the two conditions are satisfied 

1. coefficients 4>*,ip* of the system ( |2.48| ) satisfy the equalities 

$rx(s, r, x) = f3{(f)*~ 1 {s, r, x)) r , . . 

Q sx (s,r,x)=(3(r- 1 (s,r,x)) s [Z ^ } 

with arbitrary real constant /3; 
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2. there exist such four functions of one variable b(s), d(r), $o(s), @o( r ) which 
the equations 

$ r (s,r, x) = $ r (s, r, x)($ - 0), 
Q s (s,r,x) = s (s,r,x)(6 - $) 

where the functions $(s, r, x), 0(s, r, x) are defined by the formulae 

$(s, r, x) = 6(s)$ s (s, r, x) + d(r)$ r (s, r, x) + $ (s), 
0(s, r, x) = fe(s)G 5 (s, r, x) + d(r)Q r (s, r, x) + 0o(r). 

These symmetries are generated by the Lie equations 

s T = 4>*(x, t, s, r)s x + b(s), r T = ip*(x, t, s, r)r x + d(r) 

where the functions <p*,ip* are defined as follows 
if (3^0 

0*(x, t, s, r) = 0*(s, r, x)< a(s, r) exp [(3(t+ 



x 1 > 

+ / r, x)dx)} + +-^$(s, r, x) 

o p 



and if /5 = 



•0*(x, t, s, r) = ip*(s,r,x)Y:(s,r) exp[j3(t+ 
p 



0*(x, t, s, r) = 4>*(s, r, x)< a(s, r) + / 0(s, r, x)dx— 

I o 

-$(s, r) t + /0*~ 1 (s,r, x)dx 
o 

^*(x, t, s, r) = r/>*(s, r ; #n c ( s ) r ) + / r 5 x)dx— 

I o 

— G(s,r) t + /■?/;*~ 1 (s,r, x)dx 
o 



where the functions <fr, ip are defined as follows 



<j)(s,r,x) = b(s)[4>* l (s,r,x)] s ++d(r)[(f)* 1 (s,r,x)) r , 
■0(s,r, x) = b(s)[ip*^ 1 (s,r, x)] s + +d(r)[ip*~ 1 (s, r, x)] r . 



satisfy 
(2.51) 

(2.52) 
(2.53) 



(2.54) 



(2.55) 



(2.56) 



Here the functions a(s,r),c(s,r) form an arbitrary smooth solution of the linear system 
fl2.17| ) with the coefficients $ r (s, r, 0), Q s (s, r, 0) which are obtained from the formulae 
at x = 0. 



2.4 Invariant solutions and linearizing transformations for sys- 
tems with an explicit dependence on i or x 

Let the system (|2.1| ) which may explicitly depend on t satisfy the condition 1 of the 
theorem 2.1 and hence possess an infinite set of homogeneous in derivatives hydrodynamic 
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symmetries ( |2.20| ). For the corresponding invariant solutions we have s T = r T = and 
taking into account the condition s x r x ^ we obtain <f)(x,t,s,r) = 0, ip(x, t, s, r) = 
0. Using here the formulae Q2.21D , ( |2.22| ) for <j>,ij} with cq = 1 we obtain the following 
conditions for invariant solutions 
if (3^0 



a(s, r) + exp j— j3 
c(s, r) + exp < —f3 



x + J (f>(s, r, t)dt 
o 
t 

x + J ip(s, r, t)dt 
o 



0. 
0. 



(2.57) 



and if P = 



a(s, r) + x + J 4>(s, r, t)dt = 0, 
o 
t 

c(s, r) + x + J i[)(s, r, t)dt = 0. 
o 



:2.58) 



Theorem 2.4 (see |13|). Let the coefficients <f>,if) of diagonal system ( p.l|) satisfy the 



condition 1 of the theorem 2.1. Then any solution of the system (|2.57|) for (3 ^ or of the 
system ( |2.58| ) for (5 = is also a solution of the system ( |2.1| ) if the following conditions 
are satisfied 
for ^ 



t t 
In \a{s, r)\+pj <j>{s, r, t)dt) s ■ [ In \c(s, r)\+pj ip{s, r, t)dt]\ 



(2.57) 



^ 0, (2.59) 



and for /3 = 



(s,r)+ / 0(a,r,t)dt] a -[c(a,r)+ / V(s,r,t)dt] 



(2.58) 



(2.60) 



And vice versa any smooth solution s(x, t),r(x, t) of the system ( |2.1| ) can be obtained 
from the systems (|2.57|) or ( p.58|) if the condition ( p.59|) or (|2.60| ) is met in the vicinity of 
any point (xo,t ) where the condition s x (xo,t ) ■ r x (xo,to) ^ is satisfied. 

Remark 2.6 (see |i~3f). Equations ( p.57|) for p ^ and (|2.58| ) for P = determine 
in an implicit form the linearizing point transformation of the nonlinear system ( |2.1| ) 
with explicit t-dependence which satisfies the condition 1 of the theorem 2.1. Indeed, 
a search for solutions s(x, t), r(x, t) of the system ( |2.1| ) reduces to a search for solutions 



a(s,r),c(s,r) of the linear system ( |2.17| ) with variable coefficients. In particular, if the 
coefficients 0(s, r), ip{s, r) of the system (|2.1| ) are not explicitly t-dependent then the 
condition 1 of the theorem 2.1 is satisfied with p = and the equations ( |2.58|) coincide 
with the classical hodograph transformation J5| 



a(s, r) + x + t(p(s, r) = 0, c(s, r) + x + ti[)(s, r) = 



(2.61) 



Consider now the system ( |2.48| ) with explicit x-dependence. Let it satisfy the condition 
( |2.50| ) of the theorem 2.3 for possessing an infinite set of the hydrodynamic symmetries 
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( |2.53| ) with b(s) = d(r) = 0. Then for invariant solutions determined by the equations 
s T = r T = with the condition s x r x ^ we obtain 

4>*(x, t, s, r) = 0, ijj*(x,t,s,r) = 0. (2.62) 

Here the functions <f)*,ip* are determined by the formulae ( |2.54| ) for (3 ^ and ( |2.55| ) for 
P = with $ = $ (s), O = o (r), = ijj — 0. Let the existence conditions for the 
implicit vector- function determined by the equations ( |2.62| ) be met. They are similar to 
the conditions ( p39|) , fll60l) . 



Then the equations fl2.62| ) determine a linearizing transformation for the system 



with explicit dependence on x and a theorem analogous to the theorem 2.4 is obviously 
valid. 



2.5 Higher symmetries of diagonal two-component systems 

Higher symmetries of the second order are generated by the Lie equations ( |2.3| ) with 
N = 2 



T T 9\-^i tj Sj Tj S X i fx} §a 



(2.63) 



which are compatible with the system (|2.1|) . 

Definition 2.2. We call the system ( p.l|) generic system with respect to second order 
symmetries if its coefficients <f), ip do not satisfy the constraints 



$(s, r, t) = In 



c(s, t)(j) + d(s, t) 



9(s, r, t) = In 



Vv 



G(r,t)ijj + H(r,t) 



(2.64) 



$ t (s,r,t) = A(s,t)0 + J B(s,t), e t (s,r,t) = E{r,t)il> + F{r,t) (2.65) 

with arbitrary smooth functions A, 5, c, <i, £7, F, G, and the functions <3>, G defined by 
the equations (|2~8| ). 

The analysis of the symmetries Q2.63p for the generic system ( |2.1| ) produces the fol- 
lowing results. 

Theorem 2.5. A necessary existence condition for second order symmetries of the 
generic system fl2.1|) coincides with the necessary and sufficient condition for the system 
( ^.1[ ) to possess an infinite set of homogeneous in derivatives hydrodynamic symmetries, 
i.e. with the condition 1 of the theorem 2.1. 

Define the function A(s, r) by the equation 



A sr (s,r) = -$ r (s,r,t)Q s (s,r,t). 



(2.66) 



The fact that the function A has no explicit independence of t is a consequence of the 
equations (|J), ([2T2D and ( ^66]) . 

Let A(s), C(r), b(s), d(r), $o(s), ©o( r ) be arbitrary smooth functions of one variable. 
Define the functions 



$(s, r, t) = A{s){$ 2 a - $ ss + 2A SS ) + A'(s)A s + C(r)(2$ r e r + $ rr - $ f 2 ) 

+C"(r)$ r + 6(s)$ s + d(r)$ r + $ (s). 
e(a,r,t)= A>(s)e s +A(s)(2Q s <S> s + 9 SS - Q 2 S ) + C(r)(6= - 6 rr + 2A rr ) 

+C'(r)A r (s, r) + b(s)Q s + d(r)Q r + 6 (r). 



(2.67) 
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If P — in the equations ( |2.12| ), then in virtue of the condition ( |2.18j ) the functions 
$(s, r), 0(s, r) do not dependent explicitly on t. Then define the functions 4>,i/j 



(f)(s,r,t) = A(s)(2$ s 0, - SS ) + C{r)[2Q r ct) r + rr - 2$ r (0 r - + 

+C'{r)(p r + 6(s)0 s + d(r)(j> r , 
$(s,r,t) = A(s)[2^ s ^ s + *j; ss - 2Q S (^ S - <f> a )] + A'(s)^ s + 

+C{r){2Q r ip r - ij) rr ) + b(s)4> s + d(r)ip r . 



(2.68) 



Theorem 2.6. (see ||15p. A diagonal two-component generic hydro dynamic- type 
system ( |2.1| ) which may explicitly depend on time possesses an infinite set of second order 
higher symmetries with a functional arbitrariness which is locally parametrized by two 
arbitrary functions Ci(s), C2(r) of one variable iff the following two conditions are satisfied 

1. the coefficients <f),ip of the system ( |2.1| ) satisfy the equations ( [2.12| ) with an arbi- 
trary real constant (3 where the functions O are defined by equations ( |2.8|) and partial 
derivatives with respect to t are taken at constant values of s and r; 

2. there exist six functions A(s), C(r), b(s), d(r), $o(s), 6o(r) of one variable which 
satisfy the equalities 



$ r -($-e), e s = e s -(e-$) 



with the functions $(s, r, t), ©(s, r, t) defined by the formulae ( [2.671) . 
These symmetries are generated by the Lie equations 



Sr = f= A(s)°-^ + ^ 



+ 



+p h s x v(x,t,s,r) 



2A{s)$ s + b(s) 



C{r) 



C{r) 



A(s)- + C(r) 



s. 



r :r 



+ 



+ r x p(x, t, s, r) + 2C(r)© r + d(r) 



(2.69) 



(2.70) 



with the coefficients u, p defined by the formulae 
if /3^Q 



u(x, t, s, r) 



a(s, r) exp |/9 
p(x, t, s, r) = c(s, r) exp j/3 
and if (5 = 



x + J 4>(s, r, t)dt 
o 
t 

x + J ip(s, r, t)dt 
o 



u(x,t,s,r)— a(s, r) + j 4>(s, r, t)dt— 



-$(s,r) 



x + j 0(s, r, t)dt 



p(x,t,s,r) = c(s,r) + f ip(s,r,t)dt- 



-©(s,r) 



x + J ip(s, r, t)dt 
o 



(2.7i; 



(2.72) 
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where the functions <j), tp are defined by the formulae (|2.68| ). Here the integrals with respect 
to t are taken at constant values of s,r. The functions a(s , r) , c(s , r) form an arbitrary 
smooth solution of the linear system ( [2.17| ). 



Remark 2.7. For A(s) = C(r) = the higher symmetries ( 2.70 ) reduce to the hydro- 



dynamic symmetries ( |2.14| ) and the coefficients z/, p coincide with 0, ip. 
2.6 First order recursion operators 

An effective description of an infinite set of symmetries of any order is obtained by means 
of recursion operators which by their definition map any symmetry again into a symme- 
try. Here we consider recursion operators which belong to the class of matrix-differential 
operators 

R = A N D^ + A N ^ 1 D^- 1 + ... + A 1 D x + A , (2.73) 

where Ai are n x n matrices for n-component system and here n = 2. The matrices Ai 
may depend on s, r and their derivatives of a finite order with respect to x. 

Definition 2.3. If A^ ^ 0, then the integer N is called order of the recursion operator 



Here we consider the case N — 1. 
Define the functions 

S(s, r) = r, 0) = A(s)$ s (s, r, 0) + C(r)$ r (s, r, 0) + $ (s), 
T(s, r) = G(s, r, 0) = A{s)Q s {s, r, 0) + C{r)Q r {s, r, 0) + 9 (r) 

where A(s), C(r), $o( s )> ©o( r ) are arbitrary smooth functions of one variable. 



;2.74) 



Theorem 2.7 (see [15|). Let the system (|2.1|) satisfy the condition 1 of the theorem 
2.1. Then a first order recursion operator for the system ( [2.1|) exists iff there exist such 
functions A(s) and C(r) which satisfy the conditions 

S r (s, r) = $ r (s, r, 0)(S - T), T s (s, r) = Q t (s, r, 0)(T - S) (2.75) 



with the functions S(s, r), T(s, r) defined by the formulae Q2.74Q with &o(s) = 6o(r) = 0.[] 
This recursion operator is given by the following formula 



A(s)D ie [$(s, r, t)}, -$ r (s, r, t)[A(s)r s - C(r)sJ 

e,(s, r, t) L4(s)r x - C(r) S J C(r) J D ;c [0( S , r, *)] 

( l/s x \ 

' I l/r x ) ' 



(2.76) 



A first example of first order recursion operator for a class of two-component systems 
was given in the author's paper 



Theorem 2.7 generalizes Teshukov's results [26| for explicitly t-dependent systems. 
Theorem 2.8. For the system ( |2.1| ) which meets the condition 1 of the theorem 2.1 
there exists a first order recursion operator R iff there exists a recursion operator R which 

1 It is equivalent to the redefinition & s {s, r, 0), 9 r (s,r, 0). 
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R 



acts on a subspace of homogeneous in derivatives hydrodynamic symmetries of the form 
( j2.20| )- (|2.22j ) with Co = 0. An action of the operator R is defined by the appropriate 
restriction of R 

( 1 \ 

1 a(s, r) exp {(3[x + / <f>(s, r, t)dt]}s x ' 

o 
t 

c(s, r) exp {(3[x + J i/)(s, r, t)dt]}r x 
o 

/ * \ 

a(s, r) exp {(3[x + / 0(s, r, t)dt]}s x ' 

o 
t 

c(s, r) exp {P[x + J if)(s, r, t)dt]}r x 
o 



V 



/ 



(2.77) 



V 

with functions d(s,r),c(s,r) defined by the formulae 



a[s, r) 



c(s,r) 



A(s) [a s (s, r) - $ s (s, r, 0)a(s, r)] - <$> (s)a(s, r) 

- C(r)$ r (s,r,0)c(s,r), 
C(r)[c r (s,r) - Q r (s,r, 0)c(s,r)] - o (r)c(s,r) 

- A(s)6 5 (s,r,0)a(s,r). 



(2.78) 



Corollary 2.5. For any smooth solution a(s, r), c(s, r) of the linear system (p.17 ) the 
functions a(s , r) , c(s , r) defined by the formulae ( |2.78| ) are also a solution of this system 
iff the condition ( |2.75|) is met. Formulae ( p. 78j ) determine a recursion of solutions of the 



system ( |2.17| ). 

Remark 2.8. In virtue of the linearizing transformations ( 2.57 ) or ( p. 58 ) a search 
for any " nonsingular" solutions of the nonlinear system ( |2.1| ) reduces to a search for 
solutions a(s , r) , c(s , r) of the linear system ( |2.17|) with variable coefficients. However, 
the integration of the equations fl2.17| ) is also a problem. The existence of a recursion 
operator and hence of the recursion fl2.78|) which allows to multiply solutions of the linear 
equations (p. 17 ) is the important property of the system (2.1). It allows us to go over 



from linearization of the system ( 2.1|) to its integration. For this to be true its coefficients 
4>,ip have to meet the condition 1 of the theorem 2.1. and the conditions (|2.75|) . 

Remark 2.9. Linear system (|2.17|) has two trivial solutions (a, c) = (—1,-1) and 
(a(s,r), c(s,r)) = ((p(s,r,0),ip(s,r,0)) (see the equations ( |2.8|) ). The recursion formulae 
( |2.78| ) map them to nontrivial solutions d(s , r) , c(s , r) of the equations ( |2.17| ). 

In particular, the first solution is mapped to the solution (|2.74|) : a = S(s,r),c = 
T(s,r). Substituting these new solutions for a,c to the equations (|2.78|) we obtain new 



nontrivial solutions of the system ( p. 17]) . Thus, we construct two infinite series of its 
solutions. Then the linearizing transformations (|2.57| ) and (|2.58|) generate two infinite 
series of exact solutions of the equations ( |2.1| ) in explicit form. 

Corollary 2.6. The transformation (|2.78| ) maps any smooth solution (a, c) of the 
system ( |2.17| ) again into a solution (a, c) of this system iff it maps the trivial solution 
(a = 1, c = 1) into some solution of the equations ( |2.17|) . 

Corollary 2.7. The recursion ( p.47| ) for solutions of the linear system ( |2.17| ) gener- 
ated by the Lie commutator with the inhomogeneous in derivatives hydrodynamic sym- 
metry coincides with the recursion ( p.78[ ) generated by the first order recursion operator 
( j2.76 ). The existence conditions for the recursion operator and for the inhomogeneous 
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in derivatives hydro dynamic symmetry of the system ( |2.1| ) coincide. If there exist sev- 
eral such symmetries of the form ( |2.24|) , then there exist several corresponding recursion 
operators of the form ( |2.76| ) with A(s) = b(s), C(r) = d{r). 

Recursion operator R of the form (|2.76|) generates an infinite set of higher symmetries 
( |2.3D of any order N. Second order symmetries are generated by a twofold action of the 
operator R on homogeneous in derivatives hydrodynamic symmetries ( 2.20Q — ( 2.22 ) with 
c 7^ where we can put c = 1, a(s, r) = c(s, r) = 0. For instance, in the case (3 = the 
formulae ( [2.20|) and (|2.22 ) give the following form of this initial symmetry 



x + 0(s, r, t)dt 



x + ip(s, r, t)dt 



(2.79) 



2.7 Second order recursion operators and higher symmetries 



Definition 2.4. We call the system ( |2.1[ ) generic with respect to second order recursion 
operators if its coefficients do not satisfy any of the following constraints 



ci(s e 



Q S r 

e7 



c 2 (r)e 



e 



(2.80) 



with arbitrary smooth functions Ci(s), C2(r). 
Define the functions 



S{s,r) 
T(s,r) 



A(s)(<S> 2 s - $ ss + 2A SS ) + A'(s)A s + C(r)(2$ r Q r + $ r 
+C"(r)$ r + b(s)$ s + rf(r)$ r + $ (s), 
A'(s)Q s + A(s)(2& s $ s + & ss - Q 2 S ) + C(r)(9 r 2 - rr 
+C'(r)A r + b(s)e s + d(r)e r + 6 (r). 



2A ........ 



(2.81) 



with arbitrary smooth functions A(s), C(r), b(s), d(r), $o(s), @o( r ) and the functions 
$(s, r, t), Q(s, r, t), A(s, r) defined by the formulae fl2.8Q and ( |2.66| ). 

Theorem 2.9. (see Jl^]). Let the generic system (|2.1|) satisfy the condition 1 of the 
theorem 2.1. Then a second order recursion operator of the form ( |2.73| ) with N = 2 exists 
for the system ( |2.1| ) iff there exist six functions A(s), C(r), b(s), d(r), $o(s), Qo(r) of one 
variable which satisfy the conditions 



S r (s, r) = $ r (s, r, 0)(S - T), T r (s, r) = Q s (s, r, 0)(T - S). 
This recursion operator is determined by the formula 



(2.82) 



( - \ 



R = (AD 2 X + BD x + F) 



— 



[2.83) 



with 2x2 matrices A, B defined as follows 

/ A(s) 



A 



\ 









C(r) 



;2.84) 
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B 



A(s) 



si 



-Q s (s,r,t) 



2A(s) L s + $ r r ± + ^- 
\ s x s x 

A(s)^-C(r) S -^ 
A{s)^-C{r)^ 

S x fx 



b(s) 



C{r)^ + 2C(r) 6 S ^ + Q r + ^ + d(r 



f ... 



and the elements fij of the matrix F are defined by the equations 



fl2 = A(s) 



^r_x_ fSxx _ j + ^ 



-C(r)A sr ^ + {A{s)[2{$ rs - $ r $ s ) - A sr ] - b(s)$ r } r x + 



+ {C(r)(2$ r 6 r + $ rr - $2) + [C'(r) + d(r)]$ r } S(B - 



/21 = C(r) 



e. 



C(r)? 



?x r x J r x 



rt 



-A(s)A sr -t + {C(r)[2(Q sr - Q s Q r ) - A sr ] - d(r)Q s } s x + 

+ {A( s )(2e s $ s + e ss - 62) + + b( s )]e s } r x + 



A( s ) r -^-C{r) S -^ 

Sx ^x 



A(s) r -^ + C(r) S -^ 



fn + fn = Hs,r,t)s x + p{A(s) 
+<$> r (s,r 7 t) 

f. i: ■ J.,-, Q(s,r,t)r x i{('(r 
+Q s {s,r,t) 



% + 2$ s ( S ,r,t) + ^ 



+ &(*)}, 



f cr.cr. 



+ 26 r (s,r,t) + 







+ 



+ 



A ( s ) r J!L + C(r)^ 



d(r)}. 



(2.85) 



(2.86) 



(2.87) 



Here $ = $(s, r, £), 6 = G(s, r, t) and the functions $(s, r, t), 0(s, r, t) are defined by the 
formulae G2.67| ). 

Theorem 2.10. For the system ( |2.1| ) which satisfies the condition 1 of the theorem 
2.1 there exists a second order recursion operator R of the form ( [2.83 ) iff there exists the 
recursion operator R which acts on a subspace of homogeneous in derivatives hydrody- 
namic symmetries of the form ( p.20| )- (|2.22| ) with Co = 0. The action of the operator R 
is defined as an appropriate restriction of R by the formula ( |2.77| ) where the functions 
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a(s,r),c(s,r) are defined as follows 



a = A(s)a ss - [2A(s)$ s + b(s)]a s - C(r)$ r c r + 
+[A'(s)A s + A(s)($ 2 s - $ ss + 2A SS ) + b(s)$ s + $ (s)]a+ 
+[C"(r)$ r + C(r)(2$ r 6 r + $ rr - + d(r)$ r ]c, 

c = C(r)c rr - [2C(r)© r + d(r)]c r - A(s)6 s a s + 
+ [C"(r)A r + C(r)(9 f 2 - Q rr + 2A rr ) + d(r)Q r + o (r)]c+ 

+[A'(s)e, + a(s)(2& s <s> s + e ss - e 2 s ) + b(s)e.]a. 



(2.{ 



Corollary 2.8. For any smooth solution a(s , r) , c(s , r) of the linear system ( 2.17|) 
the functions a(s, r), c(s, r) defined by the formulae ( p.88|) form also a solution of this 
system iff the condition ( |2.82j ) for the functions ( 2.81| ) is satisfied. The formulae ( |2.88| ) 



determine a recursion of solutions for the system ( 2.17|) . 

The Remarks 2.8 and 2.9 and the Corollary 2.6 are still valid for the recursion (T. 
Thus, again we can construct two infinite series of solutions of the system (|2.1|) starting 
from a trivial solution of the system (|2.17|) . 



Theorem 2.11. Existence conditions ( ]2.82j ) with the notation ( |2.81| ) for the second 



order recursion of solutions of the system (|2.17|) are less restrictive than the existence 



conditions ( |2.75| ) with the notation ( |2.74| ) for the first order recursion ( |2.78| ). 



Higher symmetries of the system ( |2.1|) are generated by the action of the recursion 
operator ( |2.83| ) on the homogeneous in derivatives hydrodynamic symmetries ( 2.2U )-( 2.22 ) 
with Co 7^ 0. 



Theorem 2.12. (see ||15||). Second order symmetries ( |2.70| ) of the system ( p7T 



obtained as a general solution of determining equations coincide with the second order 
symmetries obtained by the action of the second order recursion operator ( |2.83|) on the 



hydrodynamic symmetries fl2.20|) -( pi.22|) with cq ^ 0. Existence conditions for the second 
order higher symmetries and for the second order recursion operator also coincide. 

Corollary 2.9. All the second order symmetries may be obtained by an action of 
the second order recursion operator on the homogeneous in derivatives hydrodynamic 
symmetries. 

Remark 2.10. The method of calculation of higher order recursion operators developed 
1M is much more simple than a straightforward calculation of higher symmetries from 



in 



determining equations. Thus, with the suitable extension of the corollary 2.9 we see that 
this easier way of calculation of symmetries by means of the A^th-order recursion operators 
gives all higher symmetries of the same order. In particular, a squared first order recursion 
operator does not produce a general form of second order symmetries. 



3 Generalized gas dynamics equations 
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3.1 Symmetries of one-dimensional isoentropic 
gas dynamics equations 

We consider one-dimensional gas dynamics equations for the isoentropic plane-parallel 
gas flow 

u t + uu x + a 2 (p)pp x = 0, , . 

Pt + pu x + up x = 0. 

Here u(x, t),p(x, t) are gas velocity and density at the point x at the moment t, c — pa(p) 
is the speed of sound, a(p) is an arbitrary smooth function. In practice oi{p) is determined 



by the gas state equation p = P(p), where p is a gas pressure: a(p) = (l/p)yP'(p). 
Theorem 3.1. (see ||). The system ( |3.1|) can be reduced to diagonal form ( p7T 



[s,r) 



s t = 
s + r 



{s,r)s x , r t = ^(s,r)r a 



pa{p) 



s + r 



pa{p) 



(3.2) 



by the transformation to Riemann invariants s, r 



u 



a(p)dp, r = u + a(p)dp 



PO 



(3.3) 



where po is an arbitrary fixed constant. The inverse transformation is given by the for- 
mulae 



r + s 



u 



a(p)dp 



(3.4) 



Pa 



where the last equality determines p as an implicit function p = p(r — s) for any fixed 
choice of the function a(p). 

The determining equation for symmetries of system ( |3.2|) has the form 



(ID t + A)(f,gY =(0,0) 



(3.5) 



where / 



1 
1 



T means a transposed matrix, the operator of the total derivative 
D t is calculated with the use of the equations ( |3.2| ) and the operator A has the form 

s + r f ^ „ \ 



A 



pa(p), 











+ 



s x 
r T 



1 + 



2 

2a(p) 



+ pa{p) ) 

1 -1 
-1 1 



D x + 



(3.6) 



Canonical Lie-Backlund symmetries ( |2.3| ), (|2.4| ) of iVth order for = 1,2,3 were 
obtained as general solutions of the determining equation (|3.5| ) for a generic a(p) (see 
16|.[17|]) explicitly independent of x, t 
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N = 1 : 

ft = a(s, r)s x , gt = c(s, r)r x 
where a(s,r),c(s,r) is an arbitrary smooth solution of the linear system 

a r\s,r) = ——{a- c), c 8 {s,r) = -—(a- c), 



N = 2: 



Aa 2 



4a 2 



, _ _Sxx_ _ (Aft) (sx ~r x ) 2 
si Aa 2 s x r x 

- _|_ ( s x ~r x ) 2 

r 2 4a 2 s x r x 



(3.7) 



(3i 



(3.9) 



Af = 3 : 



/a 



3a' (p) 
4a 2 s x 



3s xx 



si 



2M fi _ L) 

4a 2 \ n! / 



^3 = - 



2 \a 2 
3r 2 



3a' (p) ^ 
4a 2 r? 



fx + 



1 

2 la 2 



"x 'a 

<*'(p) 
4a 2 

1 1 
— + — 



cr / as T 



n Q I T"xS XX 

. si rt } 



K s x - r x ) 
8s x r x 



a 



cr 



(3.10) 



Here every A^th-order symmetry is presented up to addition of lower-order symmetries. 
Hence all the higher symmetries have a functional arbitrariness determined by the linear 
system ( |3.8| ) since we can add to them the terms (|3.7|). 

Analysis of the determining equation (|3.5|) for second order (A" = 2) symmetries 
produces the following special choices of a(p) which lead to extensions of the set of 
symmetries. 

1. Function a(p) satisfies the differential equation 



{cp + b)a 2 {p)/a\p) 



(3.11) 



with arbitrary constants A, b, c. The second order symmetries are determined by the for- 
mula 



(As + A) f r ±- 2 ) + (Ar + C)^ 



+ a(s,r)s x , 



g2 = (Ar + C) 



i ... 



Ap) 

Aa 2 



(Ar + C) f^-2) + (As + A)^ 



+ b 



+ b 



(3.12) 



+ c(s, r)r x 



with arbitrary constants A, A, C, b and functions a(s, r), c(s, r) satisfying the linear system 
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In particular, the equation ( 3.1 1|) is satisfied for c = by the physically interesting 
state equation of a polytropic gas 



P{p) = a 2 p\ a(p) = aV7P (7 " 3)/2 (3.13) 

where a, 7 are constants . 

2. Function a(p) satisfies the condition a'(p) = 0. This is a polytropic gas with 7 = 3. 
Then the second order symmetries 

„=.,„( t ..). 

depend on two arbitrary smooth functions ipi,ip2- 
The gas dynamics equations ( |3.2| ) have the form 



St = —ss x , r t = -rr x . (3.15) 

Their general solution 

x-st = F(s), x-rt = G(r) (3.16) 

depends on two arbitrary smooth functions F and G. Hence s T = r T = and the solution 
manifold consists solely of invariant solutions. 

3. The function a(p) satisfies the condition a'(p) = — (2/p)a(p). This is the Chaplygin 
gas [^7] with the state equation 

P(p) = P --, «(p) = 4 (^o>0) (3.17) 
P P 

with the constants Pq, a. Then the set of second order symmetries depends again on two 
arbitrary functions V'i( s ; s q),' l p2{ r -, r q) where q is the Lagrangian coordinate defined by the 
equation B 

dq = pdx — pudt. (3.18) 

The diagonal form (|3.2|) of gas dynamics equations after the transformation to the 
Lagrangian coordinates q, t becomes 

s t = as q , r t = -ar q . (3.19) 

Its general solution 

s = F(q + at), r = G{q-at) (3.20) 

depends on two arbitrary smooth functions F and G. Since x T = t T = and hence q T = 0, 
we have s T = r T = as well and the solution manifold again consists solely of invariant 
solutions. 

In the last two cases the reason for the gas dynamics equations to be integrable in 
explicit form is that the extent of arbitrariness of the set of invariant solutions and of the 
general solution manifold turns out to be the same, i.e. two arbitrary functions of one 
variable. Hence all the (nonsingular) solutions are invariant. 
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3.2 Recursion operators for gas dynamics equations 

Theorem 3.2. A first order recursion operator for symmetries of the gas dynamics 
equations (EO) is given by the formula 



R 



W x 



2a 



1 -1 
-1 1 



( - \ 



V 



So-. 







(3.21) 



/ 



It commutes with the operator of the determining equation fl3.5|) on the solution 
manifold of the equation (|3.5| ) 

[ID t + A,R] = (3.22) 

where the operator D t is calculated with the use of the gas dynamics equations. 

Corollary 3.1. Operator (ggg ) raises the order of higher symmetries by one unit 
according to the recursion formula 



R(U,9nY =(/w,ta)> (iV = 2,3,...) 



(3.23) 



and thus, generates an infinite countable set containing Lie-Backlund symmetries of any 
order. 

We note the equalities 



h 

92 



R 



R 2 x 



s 
r 



(3.24) 



Corollary 3.2. A solution of the recursion relation (|3.23|) has the form 



f N 
9n 



R 



N-l 



R N x 



(iV = 2,3,...). 



(3.25) 



Formula ( |3.21|) gives us a general form of the first order recursion operator for a generic 
function a(p). If the function a(p) satisfies the equation 



a' 



)'[A(r - s) + c-a] 



-2A 



a: 



a 



(3.26) 



with arbitrary constant A, a, c, then we obtain a more general form of a recursion operator. 
In particular, the equation ( |3.26|) is satisfied by the state equation ( |3.13|) of the polytropic 
gas if the constants satisfy the equations 



.4 



il-l)(a-c)/(4a^pt m ) (Pot^O), 
a = c (p = 0). 



(3.27) 



Theorem 3.3. For equations of the polytropic gas there exists a first order recur- 
sion operator which depends upon one essential arbitrary constant d/c (po ^ 0) or 
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a I A (po = 0) and has the form 



4a 2 



(As + a)(r x - s x ) -[(As + a)r x - (Ar + c)s a 

-[(As + a)r x - (Ar + c)s x ] (Ar + c)(r x - s x ) 

/i \ 





(3.28) 



1 



V 



/ 



where a'(p) / (4a 2 ) = (7 — 3)/[8a A /7p^ 7 1 ^ 2 ] and the constants A, a, c satisfy the equations 
& 

Now consider the action of the recursion operator ( |3.21|) on the subspace of hydrody- 
namic symmetries ( |3.7|) subject to conditions ( p.8|) 



R 



a{s,r)s x 
c(s,r)r x 



a!{s,r)s x 
ci(s,r)r x 



Here the functions a\, c\ are generated from a, c by the transformation 

ai(s,r) = a s (s,r) + [a' (p) / (Aa 2 )](a - c), 
ci(s,r) = c r (s,r) + [a'(p)/(Aa 2 )}(a - c). 



(3.29) 



(3.30) 



Corollary 3.3. The transformation ( |3.3C| ) determines a recursion for solutions of 



the linear system (|3.8|) , i.e. if a(s,r), c(s,r) is its solution, then ai(s,r), ci(s,r) is also the 
solution. 

From the equation ( |3.29|) we obtain 



R 



n ( a(s,r)s x 
c(s,r)r x 



a N (s,r)s x 
c N (s,r)r x 



(N = l,2,...) 



(3.31) 



where cin(s, r), cn(s, r) is a result of iV-fold application of the transformation (|3.30|) to 
the solution a(s, r), c(s, r) of the equations (|3.8|). Then a^(s, r), cn(s, r) is also a solution 
of the equations (|3.8j ). 



3.3 Generalized gas dynamics equations, their symmetries 
and recursion operators 

In the Lie equations for hydrodynamic symmetries ( ^.7[ ) of gas dynamics equations 

s T = a(s,r)s x , r T = c(s,r)r x (3.32) 



subject to the conditions (|3.8|) we consider the group parameter r as a new time variable. 
Then we change the notation a(s,r),c(s,r) to <p(s,r),ip(s,r). We obtain systems of the 
form (|2.1|) subject to the additional constraints 



s t = 0(s, r)s x , r t = tp(s, r)r x , 
r (s,r) =^ s (s,r) = [a'(p)/4a 2 ](0-VO 



(3.33) 
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i. e. (0, ip) is an arbitrary smooth solution of the linear system ( |3.8| ) . 

Equations ( |3.33|) appeared for the first time in the papers of the author [)17|,|l8| and 
later Olver and Nutku had called them generalized gas dynamics (GGD) equations ||. 
They also have pointed out many interesting applications of these equations in physics. 

The determining equation for symmetries of any of GGD systems ( |3.33| ) has the form 
similar to the equation (|3.5| ) 

(3.34) 



(ID t + A {M )(f,gy =(0,0) 
where the operator Au,m is defined by the formula 




s, r) 
?P(s,r) 



D x + 



a '(p) i t i \ 



(3.35) 



The functions <f>i,if>i are generated from <p,ib by the transformation ( p. 30 ) 



Ms,r) = <f> s (s,r) + [a'(p)/(4a 2 )](0-^), 
Ms,r) = </v(s, r) + [a'(p)/(4a 2 )}(<f) - tb) 



(3.36) 



and hence the functions 4>\,ty\ satisfy the equations ( |3.33| ). 

Theorem 3.4. For the GGD equations ( |3.33p all the hydrodynamic symmetries 
homogeneous in derivatives and explicitly independent of x, t are generated by the Lie 
equations fl3.32| ) with the coefficients a(s,r),c(s,r) satisfying the linear system ( |3.8| ) and 
hence coincide with the hydrodynamic symmetries of gas dynamics equations. 

Thus for any GGD system hydrodynamic symmetries are generated by the Lie equa- 
tions which belong to the same GGD hierarchy but have a different time variable r. 

Theorem 3.5. Operator R defined by the formula ( 3.21| ) is a recursion operator for 
symmetries of the whole infinite GGD hierarchy ( |3.33| ). It commutes with the operator 
of determining equation ( |3.34|) on the solution manifold of the equation (|3.34|) 



[ID t + A {(P>i , ) ,R] = 



(3.37) 



where the operator D t is calculated with the use of GGD equations. 

Theorem 3.6. For generalized gas dynamics equations ( 3.33| ) with the coefficients 
<p(s,r),ib(s,r) all the hydrodynamic symmetries are generated by the Lie equations 



s T ■ 



a (s,r s 3 



c s, n?\ 



X[x + t(p(s, r))s x + co[l + £0i (s, r)s x ), 
X[x + tij){8, r)]r x + cq[1 + tifa.(s, r)r x ] 



(3.38) 



where the functions a(s,r),c(s,r) satisfy the equations (|3.8| ) and A,co are arbitrary con- 
stants. 

Theorem 3.7. Generalized gas dynamics equations with the coefficients <p(s, r),ib(s, r) 
have an infinite series of higher symmetries of any order N = 2, 3, . . . with a functional 
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arbitrariness and an explicit t-dependence. Symmetries of the order N = 2 are generated 
by the Lie equations 



s T = ai(s,r)s x . - A[l +t(j)i(s,r)s x ] + c [t<j) 2 (s,r)s x + f 2 ], , ^ 

r T = Ci(s,r)r x - A[l + ty x (s, r)r x ] + c [tifj 2 (s,r)r x + g 2 \. 



Here the functions oi, c\ and 0i, are obtained by the transformations ( 3.3U| ) and ( 3.36 ) 



from the functions a, c and <f>,if) respectively. The functions 4> 2 ,ip 2 are obtained by 
twofold application of the transformation (|3.36|) to the functions <fi, ip and the functions 
a\(s , r) , c±(s , r) form an arbitrary smooth solution of the equations ( |3.8| ). The functions 
f 2 ,g 2 are defined by the equations ( |3.9|) . (N + l)th-order symmetries (N + 1 > 3) are 
generated by the Lie equations 




a N (s, r)s. \ _ A | ( In \ + f ( Ms, r)s x \ 1 
CN{s,r)r x I \\9nJ V VN{s,r)r x I J 



+c +t *™<f> r \°* ]} (iV = 2,3, 

V ^iv+i / \ *pN+i{s,r)r x 




(3.40) 



Here the subscript denotes A^-fold application of the transformations (|3.30| ) and ( |3.36|) , 
the functions f^, are defined by the formula (|3.25| ) with their explicit form for = 2, 3 
given by the formulae (|3.9|), ( |3 . 1 0| ) and A, c are arbitrary constants. 



3.4 Noncommutative Lie-Backlund algebra associated with 
gas dynamics equations 

Denote by Xu c -\ canonical Lie-Backlund operators of hydrodynamic symmetries generated 
by the Lie equations ( |3.32| ) subject to the condition (|3.8|). Let Xn denote canonical Lie- 
Backlund operators of the iVth-order higher symmetries ( |3.25| ) (N = 2, 3, . . .) for gas 
dynamics equations. 

Theorem 3.8. Hydrodynamic and higher symmetries of gas dynamics equations 
generate infinite-dimensional noncommutative Lie-Backlund algebra in which the hydro- 
dynamic symmetries form an infinite-dimensional commutative ideal 



Xm, Xn 



0. 



Xn, X( a ,c) 



X, 



(a N ,c N )j 







(M,N 



2,3,...). 



(3.41) 



Here the functions a(s,r),c(s } r) and a(s,r),c(s,r) satisfy the equations 
Theorem 3.9. Let N = 2,3, . . . . The generalized gas dynamics equations have a 
common infinite series of higher symmetries ( |3.40| ) of (iV+ l)th or larger order iff the right- 
hand sides of these equations differ only by a term belonging to a kernel of the operator 
R N , i.e. (pNiipN coincide for all these equations. If the right-hand sides (j)s x ,ipr x of GGD 
equations ( |3.33|) belong themselves to a kernel of the operator R N , then for all such 
equations the operators Xr^) commute with all the operators X n of higher symmetries 
( ggg ) of the order n = N,N + 1,... . 

Corollary 3.4. Commutative symmetry subalgebras for gas dynamics equations are 
generated by the operators X/ a A of those hydrodynamic symmetries whose characteristics 
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(as x , cr x ) belong to a kernel of some integer degree R N of the recursion operator (N > 2), 
and by the operators X n of higher symmetries ( |3.25| ) of the order n larger or equal to N. 




Thus, a problem of constructive description of the kernel of operator R N naturally 
arises. It is solved by means of the inverse recursion operator R~ l 



R~'= n J U + 5 _ :\D- X X -AD^ (3.42) 



where D~ x = J dx is the operator of integration with respect to x at a constant value of 
t with the integration "constant" c(t) depending upon t. Hence a kernel of the operator 
R has the form 

*-'(o )=*(<)( (;:*)■ < 3 -«) 

A kernel of the operator R 2 is given by the formula 




y ' \ {u + pa(p))r x V 2K \ [Ja 2 (p)dp + ua(p)]r x 



Here Cj(t) are arbitrary smooth functions. 



(3.44) 



3.5 Lax representation and invariant solutions of generalized gas 
dynamics equations 

Formula Q3.22j ) gives the Lax representation of the gas dynamics equations 



dR/dt = [R,A]. (3.45) 



Here the recursion operator (|3.21 ) and the "stationary part" ( |3.6|) of the operator of the 



determining equation form the Lax pair of the Ibragimov-Shabat type p9-j3l"H. Equation 
(p.37| ) also gives Lax representation for any generalized gas dynamics equations ( |3.33|) 

dR/dt = [R, A {M) ] (3.46) 



with the recursion operator ( p.21| ) as L-operator and the operator ( |3.35| ) as A-operator. 



In 1982 L.D. Faddeev and P.P. Kulish in a private communication pointed out to 
the author that these Lax pairs were degenerate, i.e. the mapping of a potential to 
scattering data was singular, so that the method of inverse scattering transform could 
not be applied. However, it turned out that it is not a deficiency of GGD equations but 
the indication to a possibility of linearizing these equations in a classical sense by the 
hodograph transformation and not by the inverse scattering transform. 

Nevertheless, the Lax representation for these equations turns out to be useful for 
construction of their invariant solutions by means of the inverse recursion operator (|3.42|) 

E3- 
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At first consider the solutions of the gas dynamics equations ( |3.1| ) which are invariant 
with respect to the higher symmetries given by the formula (|3.25|) 

C/V, 9n) T = R N ~\l, 1) T = (N = 2, 3, . . .) (3.47) 
Define the matrix U 

U=\ U l a2{p)dp I (3.48) 



p u 



and the matrix-integral operator K 



K = D^U X = D-^ U ^ a2{p) / x y (3.49) 

Define the operator K and its degrees K n with an integer n so that all the constants 
Cj of integrations D' 1 = J dx do not depend on t. 

Theorem 3.10. For any iV = 2, 3, . . . the equalities 

r)— ( /a2 i p)<ip ) +c -(p) + 

/ u 2 




u J a 2 (p)dp 
hfdpfa 2 (p)dp 



+ C27V-5 

x — ut ' 

-pt J 



-^fa 2 (p)pdp \ + = (3>50) 
V M ^ / 



determine an infinite series of exact solutions u = u(x,t),p = p(x,t) of gas dynamics 
equations ( |3.1|) which are invariant with respect to higher symmetries ( |3.25| ). They are 
given in the form of an implicit function t = t(u, p),x = x(u, p). 
For N = 1 the formula (|3.50| ) have the form 



C\ \ ( x — ut 
c 2 J [ -pt 



(3.51) 



and gives a trivial solution for which u t + uu x = 0, p t + pu x = 0, i.e. the velocity and 
density u and p are constant in the Lagrangian frame which moves together with a gas 
particle. 

For N = 2 the formula ( |3.50| ) can be taken in the form 

u-u = pit - t ), 
p(t - t ) 2 - J P'{\p)- 2 dp = x-x - u (t - t ), (3 ' 52) 

PO P 

where p = p/\ and A, w , ^o, ^o, Po are constants. It corresponds to the motion of a piston 
in a gas flow after the explosion. 

Now we consider invariant solutions of generalized gas dynamics equations expressed 
through variables u, p 

u t = w(u, p)u x + v(u, p)a 2 (p)p x 



'x ■ 



pt =v{u,p)u x + w{u,p)p x l "" h>> 
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where the coefficients v, w satisfy the linear system 



w u 



w p = a 2 (p)v u . 



(3.54) 



The corollary 3.4 implies that the higher symmetries fl3.25| ) of gas dynamics equations of 
the orders N Q , Nq + 1, ... are also higher symmetries of those GGD equations (|3.33|) whose 
right-hand sides belong to a kernel of the recursion operator R N °, i.e. 4>n = ipN — 0. We 
give now their explicit form allowing them to depend explicitly on t 



D 



a>2N -l(t) 
0>2N (t) 



Pa 



a 2 (p)p z 



u. 



+a 2 jv -2(i) 1 



j K N -l ^ 

fa 2 (p)dp 



ai(t) 
a 2 (t) 



u x a 2 (p)p x 
Px u x 

u 



it. 



+ O'2N -3(t)\ +«2Af -4(^) - 



P 



(3.55) 



/ 



u fa 2 (p)dp 



ir 



\ y + JdpJu 2 (p) d P 



+ «27V -5(^) 



/ 



V 



U 



+ Ja 2 (p)pdp 



+ . 



up 



Here in the definition (|3.49| ) of operator K the "constants" ctj(t) of the integration with 
respect to x may depend on t and are assumed to be arbitrary smooth functions. 

We consider the solutions of GGD equations (|3.55|) which are invariant with respect 
to the higher symmetries ( |3.25| ). 

Theorem 3.11. For any N = N , N + 1, . . . (N > 2) the equations 



JV-l 



K 



+C2N-4 



Cl 
C 2 

V 2 

V 2 

x 





_ / C27V-1 

u J a 2 (p)dp 



+ C2N-2 



fa 2 (p)dp 



+ 



u 
P 



+ JdpJa 2 (p)dp 

) + f ( a 2N -l(t) ) 

t I 



C2N-5 




+ C2N-3 



Ja 2 (p)pdp 
up 



u \ 
P J 

+ ■ 



+ 



dt + 



J o 



u 



(3.56) 



u J a 2 (p)dp 



u 



+ 



V 



v 



y + Jdpja 2 (p)dp J o 



/ a 2 N -i{t)dt+ 



1 1 



y + / a 2 {p)pdp 
up 



fa 2 N -5(t)dt + . . . 



determine an infinite series of exact solutions u = u(x,t),p = p(x,t) of any generalized 
gas dynamics equations of the form ( |3.55| ) (with arbitrarily fixed integer N Q ) which may 
be explicitly t-dependent. These solutions are invariant with respect to the iVth-order 
higher symmetries fl3.25p of these equations. Here the functions Oj(t) and the constants 
q with i < must be put equal to zero. The definition of the operator K is the same as 
in the formula (|3.50|) . 

For the gas dynamics equations (|3.1| ) we have a2N -s(t) = —1, a-i(t) = for i ^ 2N — 3. 
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4 Separable two-component Hamiltonian systems 

4.1 Hamiltonian structure of generalized gas dynamics 
equations 

We consider two- component hydro dynamic-type Hamiltonian systems of the form 



(4.1) 



where o\ = ( J J J is the Pauli matrix. Here H(u,p) is Hamiltonian density of the 

oo 

hydrodynamic-type which corresponds to the Hamiltonian 7i — J H(u,p)dx. For short 

— oo 

we call H(u, p) also a Hamiltonian. 

Equations ( |4.1| ) take the form of the Hamilton equations 



u 



u 

p 



H 



(4.2) 



with the Poisson bracket of the hydrodynamic-type |! 

{H,h} = (h u ,h p )D x a 1 (H u ,H p ) T . 
Define the Hamiltonian matrix 



(4.3) 



H 



Hpu H pp 
H U u H U p 



(4.4) 



Then the equations ([4.1] ) take the form 



u 
P 



H 



u 
P 



(ID t - HD X 



u 
P 








The gas dynamics equations ( ^.1[ ) have the Hamiltonian form (O) 
Hamiltonian 



H(u,p) 



p p 

.2 In i f J„ f „2l 



(4.5) 
with the 

(4.6) 



pu /2 + J dp J a (p)pdp 
o o 

The same is valid for the generalized gas dynamics equations ( |3.55| ) whose Hamiltonians 
H(u, p, t) are defined by the equality 



k n q ( o,i(t) 
{ a 2 (t) 



(4.7) 
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with the operator K defined by the formula ( |3.49| ) and they may depend explicitly upon 
t. For example the equation (|4.7|) for Nq = 2 becomes 



H(u,p,t) = ai(t) 
+a 2 (t) 



p p 

pu 2 /2 + Jdpfa 2 (p)pdp 



o o 
p p 

ii 3 /3! + u Jdp ja 2 (p)dp 
o o 
p p 

2 In , r j„r„,2 



+ 



+ 



(4- 



+a 3 (t)up + a 4 (t) 

+a 5 (t)p + a 6 (t)u. 



u 2 /2 + J dp Ja 2 (p)dp 
o o 



+ 



4.2 Separable Hamiltonian systems 

oo oo 

Definition 4.1. We say that Hamiltonians J H dx and / hdx commute if the Poisson 

— oo — oo 

bracket ( f4.3| ) of their densities is the exact derivative with respect to x of some function 

{H(u, p, t), h{u, p,t)} = D x [Q{u, p, t)}. (4.9) 

Then we also say that Hamiltonians H and h commute and the Hamiltonian matrices 
H and h defined by the equation (fO|) commute also: [H, h] = 0. 

Theorem 4.1. For any hydrodynamic-type Hamiltonian H(u,p) there exists an in- 
finite set of Hamiltonians h(u,p) which commute with it and with each other. They are 
arbitrary smooth solutions of the wave equation 



Hpph uu hppH uu 



(4.10) 
(4.11) 



or for H pp ^ 

h uu - V(u,p)h pp = 

with V(u,p) = H uu /H pp . 

Definition 4.2. Those systems (|4.5|) for which the wave equation ([4.1 1|) admits a 
separation of variables u, p, i.e. V(u,p) = j3 2 (u) / a 2 (p) , are called separable Hamiltonian 
systems [|6|,|18[|. Their Hamiltonians H satisfy the equation (H UU H PP ^ 0) 



'l/(5 2 {u))H uu = {l/a 2 {p))H i 



pi' 



and just the same for all the Hamiltonians h commuting with H 

(1/P 2 (u))h uu = {l/a 2 {p))h pp . 



(4.12) 



(4.13) 



Gas dynamics equations (3J) and GGD equations (|3.33|) with the Hamiltonians (4i)) and 
( (4.7[ ) are examples of separable systems with f3 2 {u) = 1. 

Definition 4.3. If H(u,p) is a Hamiltonian and the ratio H uu /H pp = l/a 2 (p) is a 
function of p only, then H is called Hamiltonian of generalized gas dynamics. In this case 
P(u) = l. 

Many interesting physical applications of the separable Hamiltonian systems are given 
in the paper 0. 
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4.3 Second order recursion operator and Lax representation for 
separable Hamiltonian systems 

We show that separable Hamiltonian systems obtained by the Manin's construction pos- 
sess higher symmetries, a recursion operator, the Lax representation of the Ibragimov- 



Shabat type and good integrability properties [|I8|,|32 
We introduce the notation 



d u = d/du, d p = 3/ dp, d u 



-i 



du, d 



-i 



dp 



(4.14) 



where the integration with respect to one variable is performed at a constant value of an- 
other variable. The Manin's construction gives rise to two fundamental series of mutually 



commuting Hamiltonians [18 



n=0 



#(2m-l)( 1)0 ) = £ (d-2^2( M))m -n-l (9 -2 a 2 (p)) n Mp) 



n=0 



lib 

^(0,1) = E (d- 2 f3 2 (u)) m ~ n (d; 2 a 2 (p)) n u, 



n=0 



^(2— i)( ,l) = E (d- 2 (3 2 (u)) m ~ n (dz 2 a 2 (p)y 



(4.15) 



(4.16) 



?1=0 



with m = 0, 1,2,... and ^-^(l.O) = 0, if^O, 1) = I. Here every operator fl*" 1 ,^ 1 
acts on all the factors standing to the right of it, e.g. 



u u 

du 





u u 



(3 2 (u) J du J (3 2 (u)du 





du. 



(4.17) 



For arbitrary constants c\, c 2 we define 

H^ N \ Cl ,c 2 ) = Cl H^ N \l,0) + c 2 H (N \0,l). (4.18) 
The basic Hamiltonians for the GGD equations ( |3.55| ) are obtained at (3(u) = 1 



# (A °(i,o) = £ 



T i U N ~ 2n 



o(N-2n)\ 



(d; 2 a 2 (p)Tp, 



t JV+ll 

„,A r +l-2n 

ifW(0,l)= V — - rr(9r 2 a 2 (p))". 

V ; (iV + l-2n)! v p vp;; 

We list explicitly several Hamiltonians from the series fl4.15| ), (4.16) 

#( )(l,0)=p, HW(l,0)=up, 
HV{l,Q) = {d-^\u) + d-*a\p))p, 
H^(1,0) = (d~ 2 p 2 (u) +d; 2 a 2 (p))up; 



(4.19) 



(4.20) 



(4.21) 
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tf(°)(0, 1) = u, «(0, 1) = d~ 2 p 2 (u) + d^ 2 a 2 (p), 

H^(0, 1) = (d- 2 p(u) + dfa\p))u, (4.22) 
# (3) (0, 1) = (d- 2 (3\u)r + d-^\u)dfa\p) + (d; 2 a 2 (p)) 2 . 

In particular, for the Hamiltonians Q4.1S ), ( [4. 20|) of the generalized gas dynamics we 
obtain |J 

HW{1, 0) = pu 2 /2 + Gi(p), # (3) (1, 0) = pu 3 /6 + uG^p) (4.23) 



where 

and 
where 



p 



G x {p) = J{p-a)aa 2 (a)da (4.24) 
o 

ff(D(0, 1) = n 2 /2 + G x (p), # (2) (0, 1) = n 3 /6 + uG^p) (4.25) 

p 

G^p) = Up- a)a 2 (a)da. (4.26) 
o 

Consider at first the special case of the generalized gas dynamics Hamiltonians, i.e. 
(3{u) = 1. Define the integral- matrix operator K using the matrix U defined by the 
formula ( |3.48| ) 

K = d^U x = ]dx( U * <*(P)P*\ = 

1 i / a 2 fp) \ ^ 4 ' 27 ^ 

= J o U u du + C/ P dp = + 

Here the integral with respect to x is taken at a constant value of t and the integral 
operator in the second line is the curvilinear integral not depending on the integration 
path in the (u, p) plane. 

Theorem 4.2. For f3(u) = 1 the generalized gas dynamics Hamiltonians of the series 
( [4.19 ), ( OD[) an< ^ their combinations H^ N \c\^C2) defined by the formula ( |4.18j ) satisfy 
the recursion relation 

(H^\H^f = K «), Hi N )) T (4.28) 

for iV = 0, 1, 2, . . . with the recursion operator K defined by the formula (|4.27|) . 
The solution of the equation ( [4.28 ) is given by the formula 

( H W H^f = K N (^(o)^(o)) r = K » [cii C2 f (4.29) 

where = H ( - N \c 1 ,c 2 ). 

For arbitrary functions /3(u),a(p) we define the matrices 

TT - ( U d P W ^ ) TJ - ( ^» 9 > 2 (P) ) (AW) 

Ul ~ [ P d-^ 2 (u) h U2 -[ P u J (430) 
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and the matrix-integral operators 



o V P* P W u x J 
K 2 = d~ l U 2x = J dx 



(4.31) 



(3 2 (u)u x a 2 (p)p x \ _ 

\ Px U x J 

{u f ] TT j rr j n if (3 2 (u) \ _ i / a 2 (p) 
J U 2u du + U 2p dp = d- 1 [ H ^ > 1 j+d^i 1 Q w 



(4.32) 



(0,0) 

with the operator d~ l defined as in the equation (|4.27 ). Here again in the second line we 
have curvilinear integrals which do not depend on the integration path. 
Define also the operator 

K = K X K 2 = d x x U Xx d x l U 2x . (4.33) 

Theorem 4.3. The separable Hamiltonians (|4.15|) , ( |4.16 ) and their combinations 
H( N \ci,c 2 ) satisfy the recursion relation 

(Hl N+ Q,H<? + *>) T = K(HW, HM) T (4.34) 

for N — 0, 1, 2, . . . with the recursion operator K defined by the equation ( |4.33| ) . 
This equation has the solutions 

{Hf m \H^) T = K m {Hf\H^) T = K m (c u c 2 f (4.35) 

for m = 0, 1, 2, . . . and 

{Hf m - l \H^) T = K m -\H^\H^) T = K m - l U x {c^c 2 f (4.36) 

for m = 1, 2, . . . where HW = H^ N \c u c 2 ). 

Thus, all the Hamiltonians of the Manin's series (|4.15| ), (|4.16|) are generated by the 
recursion operator ( |4.33p (see jTB[). 

Consider now the inverse recursion operator for these Hamiltonians 

K- 1 = K 2 X K^ X = (U 2x ) •/;,(/'.,) •/;, (4.37) 

with the matrices 

(tu- = g 5 1 -A rfiu ; < 4 - 38 > 

f3 2 u 2 x - a 2 p 2 x \ -px u x J 

<^)- - ( 7%: ) • (439) 

Define the second order matrix-differential operator 

L = D x K~ l d x l = D x (U 2x y l D x (U lx y l . (4.40) 
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Consider the Lie equations for Lie-Backlund symmetries of the Hamiltonian system 
~5|) with the Hamiltonian H and the Hamiltonian matrix H 



u 



f(x,t,U, p,U x , p x , . . . ,U X N \ p x N) ) \ _ f f n 



V P J T V 9(x,t,u,p,u x ,p x ,...,uV r > i pW) J \9n ) 
assuming that x T = t T = 0. These symmetries satisfy the determining equation 



(4.4i; 



(/A - D X H) ( f ) = ( I J . (4.42) 



Theorem 4.4. Any separable Hamiltonian system (|4.5 ) possesses the second or- 
der matrix-differential recursion operator of the form ( [4.40 ) which satisfies the recursion 
relation 

L(f n , 9n) T = (fn+2, 9n + 2) T (n = 2, 3, . . .)■ (4.43) 



Corollary 4.1. On the solution manifolds of the systems ( J4.5Q the following equality 
is satisfied 

[ID t - D X H, L] = 0. (4.44) 
Thus, we obtain the Lax representation for these systems 

dL/dt = [L, A] (4.45) 



with the Lax pair of the Ibragimov-Shabat type |2!|, [30| where A = —D X H and the 
operator L is defined by the equation (|4.40|) . 

Corollary 4.2. For generalized gas dynamics equations we must put (3{u) = 1, 
U\ = U 2 = U and the second order recursion operator ( |4.40| ) is equal to the squared first 
order recursion operator Lq 

L = Ll L = D x (U x )-\ (4.46) 



The operator L coincides with the recursion operator ( |3.21|) transformed from the Rie- 



mann invariants s,r to the separable variables u, p and satisfies the relation 

L (fn,9n) T = (f n +i,9n+i) T (n = 2, 3, . . .)■ (4.47) 
Theorem 4.5. For the separable Hamiltonian systems (|4.5| ) with the generic functions 



(3(u), a(p) the first order recursion operator does not exist and the second order recursion 
operator (|4.40|) is not reduced to a square of a first order recursion operator or to a product 



of different first order operators. 

4.4 Hydrodynamic symmetries and higher symmetries 

Theorem 4.6 Separable systems ( |4.5| ) with the Hamiltonian H(u,p) which do not de- 
pend explicitly on t possess an infinite set of homogeneous in derivatives hydrodynamic 
symmetries. All such symmetries with or without an explicit dependence upon x, t are 
generated by the Lie equations 
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h p (u,p) 




respectively. Here h(u, p) is an arbitrary smooth solution of the equation (|4.13 ), and 
h are Hamiltonian matrices of the form 



Corollary 4.3. All the hydrodynamic symmetries Q4.49Q without an explicit de- 
pendence on x, t are separable Hamiltonian systems with the Hamiltonians h(u, p) which 
mutually commute and also commute with the Hamiltonian H(u,p) of the system (|4.5|) . 

Thus, any separable Hamiltonian system ( |4.5|) is included into the infinite hierarchy 
of commuting separable Hamiltonian flows. 

Theorem 4.7. Hydrodynamic symmetries ( 4.49Q of the system ( |4.5|) without an 



explicit dependence on x, t form an invariant subspace for the second order recursion op- 
erator L. This operator generates the recursion relation for the hydrodynamic symmetries 

= a lD J W U > P l)=U U ) (4.50) 




h lp (u,p) 

preserving the Hamiltonian structure of the Lie equations (|4.49|) and generates the recur- 



sion of their Hamiltonians 



hi{u,p) = -Jr^h uu {u, p) = -^--h pp (u,p). (4.51) 
Corollary 4.4. The hydrodynamic symmetries ( f4.49|) of the separable Hamiltonian 



system (|4.5| ) which do not depend explicitly on x, t form an infinite-dimensional commu- 
tative Lie-Backlund algebra. Its elements have a functional arbitrariness, i.e. depend on 
an arbitrary smooth solution h(u,p) of the wave equation ( [4.13 ). 



Denote by h m (u } p) the result of m-fold application of the transformation ( 4.51 ) to 
h(u,p) 

h m (u, p) = (r 2 (u)d 2 /du 2 ) m h(u, p) = (a- 2 (p)d 2 /dp 2 ) m h(u, p). (4.52) 
Then the formula ( |4.50| ) gives the result 



L m h ) =hm (4.53) 




Theorem 4.8. For any separable Hamiltonian system ( |4.5| ) all its higher symmetries 
of the even order n = 2m which are generated by the action of the recursion operators 
L m upon the hydrodynamic symmetries ([4.4<j|) are given by the formula 



+ (tH m + h m )[ n (4.54) 






for m = 1, 2, . . . . Here H m (u, p) is obtained by the transformation ( |4.52| ) out of H(u, p). 
For m = l this formula gives all second order symmetries of the equations ( |4.5| ). 
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Theorem 4.9. The separable system (|4.5|) with the Hamiltonian H(u,p) possesses 
an infinite series of higher symmetries of even orders 2m without the explicit dependence 
on t, x if for some integer N the Hamiltonian H meets the condition 



L N H 



u 
P 




H 



N 







(4.55) 



i. e. the vector H(u x , p x ) T belongs to the kernel of the operator L , and then m is larger 
or equal to N. These symmetries are determined by the formula (|4.54|) with H m = 
and all of them are common symmetries of all the systems ( |4.5|) with the right-hand side 
belonging to the kernel of the operator L N . 

For h m = we obtain the special form of higher symmetries subject to the condition 
(09) with m>N 




L m x 




f 2m (u,p,u x ,p x ,...,u x 2m \p x 2 ^) 

92m{u,P,Ux,Px, ■•-,U, 



(2m) (2m) 



(4.56) 



For (3(u) = 1 we have L — Lq where Lq is the first order recursion operator ([4.46 ). 
Then the equations (|4.5|) form the GGD system of the section 4.2 and have the higher 
symmetries of any order n > 2 larger or equal to N 






if the condition 



u 



P 




(4.57) 



(4.56 



is met. 

For the generic function (3{u) we shall denote by X h and X 2rn the canonical Lie- 
Backlund operators of the hydrodynamic symmetries (|4.49 ) and of the higher symmetries 
( |4.56 ) respectively. 

Theorem 4.10. Hydrodynamic and higher symmetries of the separable Hamiltonian 
systems (|4.5|) which do not depend explicitly on x and t form the infinite-dimensional 
noncommutative Lie-Backlund algebra 



Xfi, X~ h — 0, l^2m, SVh 

^2m, = 



Xo m . Xh 



!L r , 



(4.59) 



in which the hydrodynamic symmetries Xh form an infinite-dimensional commutative 
ideal. 

Corollary 4.5. Higher and hydrodynamic symmetries of the theorem 4.10 commute 
iff h m = where h(u, p) are the Hamiltonians of the hydrodynamic symmetries Q4.49D - 
The theorem 4.9 poses a problem of obtaining the kernel of the operator L N to solve 



the equation (|4.55 ). In order to obtain its general solution we have to allow an explicit 
dependence of the Hamiltonians on t and use the Manin's series ( 4.15Q , ( |4.16 ) of Hamil- 
tonians 

#(iV)(! 

— v, v) with v — 0, 1. 
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We denote 



N 



(t)#( 2fc )(0,l) + 

+a, {N J^l(t)H^(l, 0) + a 4 ( J v_ fc ) +4 (t)^ (2fc - 1) (0, 1)] (460) 
(iV=l,2,...) 



with arbitrary smooth functions ctj(t). 

Consider the Hamiltonian system 
explicitly dependent on t 



with the Hamiltonians (N = 1, 2, 



it 



H [ u N] (u,P,t) 



= m(u, P ,t)^ . (4.61) 



Theorem 4.11. The kernel of the operator L N coincides with the right-hand side of 
the Hamiltonian system ( f4.61|) 



-N 








m(u, P ,t) 



p 



where 



L- 1 = \J lx D- x U 2x D-\ 



(4.62) 



(4.63) 



Corollary 4.6. Let N be some integer. Any separable Hamiltonian system which has 
higher symmetries ( [4.54} ) with H m = (with no explicit dependence on t, x) of the order 
2m > 2N has the form ( [4.61| ). In particular, such a system has the higher symmetries 
( |l.5(j| ) which meet the condition m > N. 

Theorem 4.12. Any separable Hamiltonian system ( [4.5[ ) with the explicitly t- 
dependent Hamiltonian H(u, p, t) possesses an infinite set of the explicitly dependent 
on t, x hydrodynamic symmetries 



u 
P 



xl + H(u, p,t)dt + h(u, p) 



u 
P 



(4.64) 



with the functional arbitrariness determined by an arbitrary smooth solution h(u, p) of 
the wave equation ( [4.13 ). Infinite series of higher symmetries of any even order 2m for 
such a system is generated by the action of recursion operators L m on the symmetries 
(PI) 



u 
P 



L m x 



u 
P 



+ 



H m (u, p,t)dt + h(u,p) 



U: 

P 



(4.65) 



and have the same extent of arbitrariness as the symmetries (f4.64j ) . In the formulae (|4.64|) 
and ( |4.65| ) the integration with respect to t is performed at constant values of u, p. 
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4.5 Invariant solutions and linearization 



We are looking for series of exact solutions of any separable Hamiltonian system (|4.61|) that 
are invariant with respect to higher symmetries (|4. t55[ ) which have no explicit dependence 
on t, x due to the condition H^(u, p, t) — for m> N. 
Invariance conditions for these solutions have the form 




+ h(u,p) 





for m — N, N + 1, . . . with h(u,p) satisfying the equation (|4.1c 



(4.66) 

For obtaining these 



solutions we suggest the method which uses the Lax representation (|4.44| ) and the inverse 
recursion operator ( |4.63| ) fTT|,|T8l|. 

Define the Hamiltonian H^(u,p) by the formula ( f4.60|) with N = m and arbitrary 
constants q substituted for functions di{t) 



HW(u,p) = £ [c 4 ( m - fc ) + i# (2fc) (l,0) + c 4(m _ fc)+2 #( 2fc )(0, 1) + 



k=0 



+C 4(m - fc)+ 3^ (2 ^ 1) (l, 0) + C 4(m - jfc)+ 4i? (2 ^ 1) (0 ) 1)]. 



(4.67) 



Let H(u,p) denote an arbitrary solution of the wave equation ( |4.12| ). 
Theorem 4.13. Invariant solutions of the system ( |4.61|) with the Hamiltonian 
H^ N \u, p,t) which meet the invariance condition (|4.66| ) are given by the formula 



L P u(u,p) 

H uu (u,p) 




HW(u,p,t) 
H [ u N u ] (u,p,t) 



dt 



(4.68) 



if the existence conditions are met for the implicit function (u(x,t), p(x,t)) determined 
by the equation ( [4.681 ). 

This equation defines the linearizing transformation which reduces the solution of the 
nonlinear separable Hamiltonian system ( [4.6 1| ), which is explicitly t-dependent, to the 
solution of the linear wave equation ( |4.12| ) for H{u,p). 

Corollary 4.7. The formula 



m 




Hff(u,p,t) 



dt 



(4.69) 



for m = N, N + 1, . . . gives the infinite series of exact invariant solutions of the equations 
( fL61 ) subject to the condition Q4.66|) with h(u,p) = 0, i. e. invariant with respect to the 
higher symmetries ( |4.56| ). To obtain an explicit form of these solutions, the expression 
QL671 ) for H^ m \u,p) must be used. 

In the formulae ( [L68j ) and ( [4.69 ) integrations with respect to t are performed at 
constant values of u and p. 

Corollary 4.8. For the system ( |4.6 1| ) with the Hamiltonian H^(u,p) without an 
explicit dependence on t the formula (|4.68|) gives the hodograph transformation 



H pu (u,p) 
H uu {u,p) 




H [ u N J(u,p) 



(4.70) 
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which interchanges the roles of independent and dependent variables (x,t) and (u,p). 

Remark J^.l. Explicitly t-dependent Hamiltonian system (|4.61| ) is not linearizable by 
the common hodograph transformation ( |4.70| ). The formula ( |4.68| ) linearizes all systems 
of the form 



and presents a generalization of the hodograph transformation for 
explicitly time-dependent systems. 

Theorem 4.14. Consider any separable system of the form 



(|4.5| ) with explicitly t- 
dependent Hamiltonian H(u,p,t) satisfying the equation ( |4.12| ). Its invariant solutions 



with respect to the hydrodynamic symmetries fl4.64j) are given by the equalities 



x + j H pu (u, p, t)dt = h pu (u, p), 



/ H uu {u, p, t)dt = h uu (u, p) 



I H pp (u,p,t)dt = h pp (u,p) 



(4.71) 



where h(u, p) is an arbitrary smooth solution of the linear equation (|4.13|) . These formu- 
lae reduce the solution of any separable time-dependent system (|4.5| ) to the solution of 
the linear equation (|4.13|) for h(u,p). This is the linearizing transformation which gener- 
alizes the hodograph transformation for explicitly time-dependent separable Hamiltonian 

systems 



5 Semi-Hamiltonian equations 

5.1 Geometry of semi-Hamiltonian systems and 
hydrodynamic symmetries 

Consider a system of quasi-linear first order equations homogeneous in derivatives 

ul = Vi(u)ui (i = 1,2,..., n) (5.1) 

with the diagonal nxn matrix V(u) = diag(uj(w)). In the rest of the article no summation 
on repeated indices will be assumed. Here u = (u 1 , u 2 , . . . , u n ) is n-vector, u % are the 
Riemann invariants |J. 

Lie-Backlund symmetries of the iVth order of the system (|5.1| ) are generated by the 
Lie equations 

u l T = r)i(x,t,u,u x ,...,ui N) ) (i = 1,2, . . . ,n) (5.2) 

where r is the group parameter and we assume that x T = t T = 0. Denote r\ = (rjx, . . . , rj n ) T . 

Consider the operators D x , D t of the total derivatives with respect to x, t and for the 
calculation of D t use the system ( |5.1| ) 

^i + ^(<l + ^^), (5.3) 



d A ( , s * d ^ „, r , , d 
where d/dx and d/dt are calculated at constant values of u,u x k \ 
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The compatibility conditions u % Tt = u\ T of the systems (|5.1j) and ( |5.2| ) take the form of 
the determining equation for the symmetry characteristic 7] 



ID t - VD X - U x 



'd\r 

dU 



| [v] = 0. (5.5) 
v iu j) is the Jacobian matrix of the 



Here / is the unit matrix, U x = diag(u^), dV/dU 
mapping v — v(u). 

Hydrodynamic symmetries are obtained if we choose N = 1 in the Lie equations ( |5.2| 
Here we consider only hydrodynamic symmetries linear homogeneous in derivatives 



u l T = Wi(u, x, t)ul (i = 1,2, . . . ,n). 

Impose the condition 
Define the symmetrical connection coefficients associated with the system ( |5.1| ) Jl 



T 3 - 



(<!■■ ■!Jjj) v 'i : j 



(i^j^k^i). 



This connection is compatible with the non-degenerate diagonal metric 



9ii{u) 



H M 



9ij = (i^j), det(^)^0. 



(5.6) 
(5.7) 

(5.8) 
(5.9) 

(5.10) 



Here Hi(u) are Lame coefficients. 

The connection is determined by the metric 



(In y/g~^) u j = (In if, 



(5.11) 



Integrability conditions of the system ( |5.11| ) T\. uk = T l ik uj with the use of the expres- 



sions ( |5.8j ) take the form of the Tsarev's conditions 

[Vi,ui/(vj - Vi)] u h = [v itU k/ (v k - Vi)] u j (i^j^k^i). (5.12) 

Definition 5.1 (see |4]]). A diagonal system ( |5.1| ) is called semi-Hamiltonian if its 
coefficients Vi(u) meet the conditions (|5.7| ) and (|5.12|) . 

Theorem 5.1 (see 0). The diagonal system (5JJ which meets the condition (|5.7|) 
possesses an infinite set of the hydrodynamic symmetries without explicit dependence on 

< = Wi(u)ul (t = l,2,...,n) (5.13) 

with a functional arbitrariness iff the system ( p.l[ ) is semi-Hamiltonian, i. e. the condition 
( |5.12j ) is satisfied. These symmetries are generated by the Lie equations ( |5.13|) with the 
coefficients Wi(u) which form an arbitrary smooth solution of the linear system 



w 



Wi 



(5.14) 
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The coefficients T\^{u) are denned by the formula ( |5.8|) . All these symmetries mutually 
commute. The set of these symmetries is locally parametrized by n arbitrary functions 
cAu 1 ) of one variable. 



Theorem 5.2 (see [12 1). Any semi-Hamiltonian system of the form (|5.1| ) possesses 
an infinite set of explicitly dependent on x, t hydrodynamic symmetries with the same 
functional arbitrariness as in the theorem 5.1. These symmetries are generated by the Lie 
equations 

u\ = [wi(u) + c(x + tVi(u))]ui (i = 1,2, ...,n) (5.15) 

with the coefficients Wi(u) which form an arbitrary smooth solution of the linear system 
( |5.14| ). Here c is an arbitrary constant. 

Theorem 5.3 (see [|]]). Any semi-Hamiltonian system of the form ( |5.1| ) with the 
non-degenerate metric ( [5.1 0| ) is a Hamiltonian system iff the following components of the 
Riemann curvature tensor vanish: R % -- { = (i ^ j). Then the curvature tensor vanishes 
identically and the variables u % form an orthogonal curvilinear coordinate system in a flat 
(pseudo-Euclidean) space. 

Corollary 5.1. For the semi-Hamiltonian system (|5.1|) the following components of 
the curvature tensor vanish 

R) kl = 0(i^j^k^l), Ri kj = Ti j>uk -Ti ktUj = 0, (5.16) 

Remark 5.1. The equalities ( |5.17| ) give compatibility conditions for the linear system 
( |5.14| ) and they are equivalent to the Tsarev's conditions ( p,12| ). 



5.2 Invariant solutions and linearization of semi-Hamiltonian 

systems 

Consider explicitly dependent on x, t hydrodynamic symmetries (|5.15|) with c = — 1. 
Invariance conditions for solutions of the semi-Hamiltonian system ( |5.1| ) with respect to 
these symmetries with u x ^ have the form 

Wi{u) = tvi{u) + x (i — 1, 2, . . . , n). (5.18) 



Theorem 5.4 (see [|J). Let the functions Wi(u) in the equations (|5.18|) form an 



arbitrary smooth solution of the linear system ( |5.14| ). Then any smooth solution u l (x,t) 
of the system ( |5.18|) is a solution of the semi-Hamiltonian system (|5.1| ). Vice versa, any 



solution u l (x, t) of the system ( p.l| ) may be locally represented as a solution of the system 
( [5.18j ) in the vicinity of such a point (x ,t ) where the condition u z x (x ,t ) ^ is met for 
each value of i = 1, 2, . . . , n. 

Remark 5.2. The equalities (|5.18| ) determine a linearizing transformation for the semi- 



Hamiltonian system ( |5.1| ), i.e. the solution of the nonlinear system ( |5.1D is reduced to 
the solution of the linear system ( |5.14| ) for the functions Wi(u). It generalizes the classical 
hodograph transformation for the case of mult i- component systems and hence is called 
the generalized hodograph transformation H. We have shown here the group-theoretical 
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origin of the linearizing transformation: every nonsingular solution of the system (|5.1|) 
is its invariant solution with respect to the hydro dynamic symmetries ( |5.15|) and the 
extent of arbitrariness of the set of these symmetries and of the set of the corresponding 
invariant solutions is the same as for the general solution manifold of the system (|5.1|) . 
i.e. n arbitrary functions Cj(w z ) of one variable. 

Hence to obtain explicit formulae for the invariant solutions of the system ( |5.1|) one 
must solve the linear system ( |5.14 ) with variable coefficients. 



5.3 First order recursion operators 

By definition, a recursion operator R maps any symmetry of the semi-Hamiltonian system 
( p.lp again into its symmetry, i.e. any solution 77 = (771, . . . , r] n ) T of the determining 
equation ( |5.5|) is mapped again into the solution R[rj\. For this to be true it is sufficient 
for the operator R to commute with the operator of the determining equation 



ID t -VD x -U x [ — 



R 







(5.19) 



on solution manifolds of the equations ( |5.1|) and (|5.5|) . If N is any integer, then iVth-order 
recursion operator has the form ( j2.73j) . 

For the systems (|5.1|) it is convenient to search for recursion operators in a slightly 
different form. 

In particular, for the first order recursion operators we assume 



R = {AD X + B)U~ l 

where A = A(u), B = B(u, u x ) are n x n matrices. 
Define the functions (i — 1, 2, . . . , n) 



(5.20) 



5i(«) = Erjfc(«M«*) + 



(5.21) 



k=l 



which depend on 2n functions Ci(u l ),di(u l ) of one variable. 

Theorem 5.5 (see |26||). For the semi-Hamiltonian system ( |5.1| ) there exists a first 
order recursion operator R of the form ( |5.20| ) iff there exist 2n functions Ci(u l ),di(u l ) of 
one variable which meet the conditions 



(i + j) 



(5.22) 



with the functions Si(u) defined by the formula (|5.21 ). Matrix elements of this recursion 
operator are given by the formula 



R 



13 



11 

5 ijCj ( u i)(D x + 52ri k (u)ui) 



k=l 



+ ri j (u){ Cj (u j )ui-C i (u i )ui) 



— ) +d j (w>)5 i j 



(5.23) 



where 8u is the Kroneker symbol. 
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Theorem 5.6. For the semi-Hamiltonian system Q5.1Q homogeneous in derivatives 
hydrodynamic symmetries (|5.13|) with no explicit dependence on x, t subject to the con- 
dition ( [5.14| ) form the invariant subspace for the recursion operator ( |5.23| ). A restriction 
of the operator R to this subspace is determined as follows 



R 



( Wi(u)v} x 
V w n {u)u n x 



\ w n [u)u'^ 



(5.24) 



where the functions Wi(u) are determined by the formula 

Wiiu) = Ci(u l )wi^(u) + di{u % )wi{u) + 

n 

+ E T\ k {u)c k {u k )w k {u). 

k=l 



(5.25) 



Corollary 5.2. For any solution {wi(u)} of the linear system ( |5.14j ) the functions 
Wi(u) form also a solution of this system, i.e. the formula ( 5.25|) is a recursion for solutions 
of the system Q5.14D , iff the conditions ( |5.22p are met. 

The first order recursion operator ( |5.23D for the mult i- component system ( |5.1| ) was 
constructed at first by a straightforward solution of the equation (|5.19|) p6| . A new, more 
simple method for constructing the recursion operators, based on the study of symme- 
tries of the set of hydrodynamic symmetries of the system (|5.1| ), was developed by the 



author |TJ]]. Thus, the group-theoretical origin of recursion operators was discovered. A 
simple geometrical meaning of the existence conditions (|5.22|) for the first order recursion 
operators was also clarified in this paper. 

To present the last statement explicitly, we consider another orthogonal curvilinear 
coordinates {r*}, which will be specified later, with the Lame coefficients 



Hi(r) 



9ii{r) 



(5.26) 



Consider the rotation coefficients /3ji(r) of this coordinate system defined by the equa- 
tions 0,0 

//„•; I;,,";, i' /./)• (5.27) 

Let p = {(r k — r 1 )} denote the set of coordinate differences. 

Theorem 5.7 (see ||14|| ). First order recursion operator for the semi-Hamiltonian sys- 
tem ( |5.1| ) exists iff the rotation coefficients (3ji{r) of some curvilinear orthogonal coordinate 
system r l depend only on the coordinate differences p: = Pji(p)- 



5.4 Second order recursion operators 

For second order recursion operators we assume the form 

R = (AD 2 X + BD X + F)U~ 1 (5.28) 

where A = A(u, u x ), B = B(u, u x , u xx ), F = F(u, u x , u xx , u xxx ) are n x n matrices. 

Define the "connection potential" V(u) by a completely integrable (in the Frobenius 
sense) system 

\'„,r(") r;,l'j, (/ / ./!• (5.29) 
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Its solution V(u) depends upon n arbitrary functions of one variable. The integrabil- 
ity conditions for the system ( |5.29|) are met as a consequence of the semi-Hamiltonian 
property ( |5.12j ). 

Define the functions 



+ 



h k {u)= f k (u k ) [rj^itfc - v\ k ) - v\ kjU , 

+ [c k (u k ) - f k (u k m k (i^k), 
b u (u)= f^) + (Ti,) 2 - 2V uiui (u) 

-ftwvviW+ciWn+diiu*), 

n 

B i( U ) = b ik(u) 
k=l 



(5.30) 

(5.31) 
(5.32) 



which depend upon 3n functions fi{u l ) , Ci(u l ) , di(u l ) of one variable. 

Let Rik, Aik, Bik, Fik denote matrix elements of the operator R and the matrices 
A,B,F respectively. 

Definition 5.2. Let Hi(u) and /%(«) denote Lame coefficients and rotation coeffi- 
cients of the curvilinear orthogonal system {u 1 } and Gi(u l ) denote a function of u in- 
dependent of u l . Semi-Hamiltonian system ( |5.1|) is called a generic system (with respect 
to the second order recursion operators) if none of the following special cases is met for 
i,3,k = l,2,...,n 

Vi,«*(«) = 0, (5.33) 
Vi , u i(u) = (Fi^GiW/ga) e v ^, (5.34) 
\Pji/(vj - Vi)}[((3ki/Hi) ul /vi^ ul = 

= \Pki/(Vk - Vi^iiPji/HA^/vi^i (5.35) 
Theorem 5.8 (see (TJ). 

For the semi-Hamiltonian generic system (|5.1|) a second order recursion operator R of 
the form ( |5.28| ) exists iff there exist 3n functions fi{u l ), Ci(u l ), di(u l ) of one variable which 
meet the conditions 

lh,r[«) r>;) (i^j) (5.36) 

with the functions B^u) defined by the formulae ( |5.30| )-( |5.32| ). Its matrix elements are 
given by the formulae 



R lk = (A ik D 2 x + B ik D x + F ik ){l/u k a 

A ik = <WiK)A4 
B ik = n k [fk(u k )(ui/u k x ) - j\(u l )(u k x /u x )] (i t k), 



Ba 



-fi(u>lJ[(K) 2 } + 2[/ 4 (t**)/t*a E nM+ 



+2/ i («*)rj i + £*(«*), 



(5.37) 

(5.38) 
(5.39) 

(5.40) 
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F lk = fiMPM/v*,) [(uU<) - (u k xx /u k )] + 

+Hu%{u x f/ui\Y\ k Y k kl - fi{u%{ulf/^ x ]Y\ Kuk + (ry 2 ]+ 
+{/ i (^)(r l yt-2r 



ll.W 

k 



[2/ 1 (« ! )r« ! + Ci («f: fc K+ 

k i -p* "pi 



+b lk (u)u x - fi(u*)(u k ju x ) ^ui{v\ k v% + r;/. ; , 
+< E[fk(u k )(ui/u k x )r ik r k kj - ' f-ivP^/v^ryr^ 



1 + 

(i ? k), 



F u = Bi(u)u x - F ik- 

k+i 



(5.41) 



(5.42) 



Theorem 5.9 (see Jl4]]). For the semi-Hamiltonian system ( |5.1| ) the homogeneous in 
derivatives hydro dynamic symmetries (|5.13|) with no explicit dependence on x, t, subject 
to the condition (|5. 14 ) , form the invariant subspace for the second order recursion operator 
( |5.37| ). Its action on this subspace is determined by the same formula (|5.24j ) from the 
section 5.3 but with the different definition of functions 



+ E f k {u k )Y\ h w ku k + E hk{u)w k 

k^i k=l 



(5.43) 



where the functions b ik {u) are defined by the formulae (|5.30|) , (|5.31|) . 

Corollary 5.3. For any solution {wi(u)} of the linear system ( |5.14| ) the formula 
( p.43| ) gives also a solution {wi(u)} of this system iff the conditions ( |5.36|) are met. Thus, 
the formula ( p.43[ ) is a second order recursion for solutions of the system ( p,14[ ). 

Theorem 5.10 (see ]14]]). If the first order recursion operator exists then there 
also exists the second order recursion operator equal to the squared first order recursion 
operator. The inverse is not true, i.e. the existence conditions ( |5.22j ) for the first order 
recursion operator do not follow from the existence conditions (|5.36|) for the second order 
recursion operators. 

This means that the existence conditions for the second order recursion operator are 
less restrictive than for the first order operator. 



5.5 Generation of infinite series of exact solutions 

To obtain explicit formulae for invariant solutions of the system ( p.l\ ) one must search for 
solutions of the linear system ( |5.14| ) and substitute these solutions for the set of functions 



{wi(u)} in the linearizing transformation ( |5.18D . Existence of a recursion operator for 



the semi-Hamiltonian system ( |5.1| ) is the additional constraint which makes it possible to 
obtain particular solutions of the system (|5.14|) . 



The linear system ( |5.14| ) has two trivial solutions 

Wi = l, Wi = Vi(u) (i = 1,2, ...,n). (5.44) 

They serve as initial elements for the generation of infinite series of nontrivial solutions 
by recursion operators. 

In particular, assume that the first order recursion operator (|5.23|) exists for the system 
(|5.1|), i.e. the conditions (|5.22|) are met. It generates the recursion ( |5.25|) for solutions of 
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the system ( |5.14j ) 

n 

Wi{u) = (-RiH)i = ^2(Ri)ik{w k 
with the first order recursion operator 

d 



fe=i 



(Ri)ik — 8 



ik 



Ciiu^—j + di(u l 
ou 



+ Y\ k (u c fc (u K 



(5.45) 



(5.46) 



Two trivial solutions ( |5.44|) are mapped by the operator Ri to the nontrivial solutions 
of the system ( |5.14| ) 



WAU 



(J2i[l])< = Y. T Uu)c k (u k ) + d i (u i ) = Si( U ) 



(5.47) 



k=l 



(5.48) 



Wi(u) = (Ri[v))i = Ci(u l )v itUl (u) + di(u % )vi(u) + 

n 

+ E T\ k {u)c k {u k )v k {u). 

k=l 

Substituting these expressions for Wi{u) to the equations ( |5.18| ) we obtain the explicit 

(5.49) 



formulae for nontrivial solutions of the system (|57T 

n 

J2 T\ k {u)c k {u k ) + diitf) = tvi(u) + x 



k=l 



Ci(u l )v i:Ul (u) + di(u l )vi(u) + £ T\ k (u)c k (u k )v k (u) 

k=l 

= tviiu) + x {i = 1, 2, . . . , n). 



(5.50) 



These equalities determine the exact solutions u l = u l (x } t) of the system ( [5.1D as implicit 
functions. 

The action by powers of the operator Ri on the trivial solutions (|5.44j) generates the 
explicit formulae for two infinite series of exact invariant solutions with N = 1,2,... 

(R?[l])i = tVi(u) + x, (5.51) 

(R?[v})i = tVi(u)+x. (5.52) 

Assume now that the less restrictive existence conditions (|5.36| ) for the second order 
recursion operator are met. Then there exists the recursion ( 5.43 ) for solutions of the 

(5.53) 



system (|5.14j ) 



Wi(u) = (R 2 [w])i = ^2{R 2 )ik[wk 
with the second order recursion operator 



k=l 



fi{ul) W? + [ft{ul)rM + Ci{ul)] 

+f k (u k )T\ k {u)-^ + b lk (u). 



+ 



(5.54) 
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Here the functions b ik (u) are defined by the formulae ( |5.30| ), (|5.31|) . 

The trivial solutions (|5.44 ) are mapped by the operator R 2 to the nontrivial solutions 
of the system ( |5.14j ) and via the equations ( |5.18| ) to the corresponding solutions of the 
system (|5.1| ) 



(R 2 [l])i = E b ik (u) = Bi(u) = tvi(u) + x, 

k=l 

(R 2 [v])i = Mu^Vi^iu) + ftfi^Ti^u) + ciiu 1 )} v i>u i(u)+ 
+ E fk{u k )T\ k (u)v k)U k(u) + E b ik (u)v k (u) = tvi(u) + x. 

k^i k=l 



(5.55) 



The powers R 2 generate two infinite series of exact invariant solutions out of initial 
solutions ( ggg ) with N = 1, 2, . . . 

(i%[l])i = tVi{u)+x, (5.56) 

(R%[v])i = tVi(u)+x. (5.57) 

We can also use linear combinations of solutions of the two series for the functions Wi(u) 
in the equations ( |5.18|) 

ci(^[l])i + c 2 (R™ 2 [v])i = tVi(u) + x (5.58) 

with any integers N, M, the operator Ri or R 2 and arbitrary constants Ci, c 2 . 

Further generalization is obtained if we substitute for the characteristic rj d = x + tvi(u) 
of the dilatation symmetry group in the formulae ( |5.15| ), ( |5.18| ) the results of the action 
on r) d of the operators i?f or R\ where L is any integer. Then we obtain the formula 

ci(^[l])i + c 2 {R™ 2 [v])i = t{Rx,M)i + ^(^ 2 [l])t- (5-59) 

In the cases when L > N and (or) L > M this formula is equivalent to using negative 
powers of the operators Ri or R 2 in the equations ( |5.58|) . 

Remark 5.3. The operators R\ and R 2 coincide with the first and second order sym- 
metry operators for the linear system ( |5.14| ) essential for the separation of variables in 
these equations (see |34],[35|]). Solution of the system (|5.14| ) by means of the separation of 
variables would mean solving completely the original nonlinear system ( p.l[ ). Therefore 
the linearizing transformation ( |5.18j ) may be considered as the extension of the method 
of separation of variables to nonlinear systems "rich in symmetries". 



5.6 Higher symmetries of semi-Hamiltonian systems 

Theorem 5.11 (see |12[). All second order symmetries of the semi-Hamiltonian system 



(5.1) are generated by the second order recursion operator ( 5.37 ) out of the hydrodynamic 
symmetries ( |5.15| ) (with c = 1). The corresponding Lie equations have the form 



u\ = T} 2>i = E Rij[u J x (x + tVj(u) + Wj(u))] = 

= -/ 4 (« i )«C/[(4) a ] + (fi( ui )/ u i)p.~ r ij u i+ 

+< E Vy^/ui + 2/^)4 + c t {/) + ui[xBi(u)+ 
+t(R 2 [v})i + Wi(u)} (i = l,2,...,n) 



(5.60) 
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where the set of functions {wi(u)} is an arbitrary solution of the linear system ( |5.14j ) 



The existence conditions ( |5.36|) for the second order symmetries and for the second order 
recursion operators coincide and must be met by a choice of functions fi(u l ), Ci(u l ), d^u 1 ). 
In the case when B^u) ^ and (_R 2 M)j 7^ these symmetries explicitly depend on x, t. 

The action of the powers R N of the recursion operator with N — 1, 2, . . . on the same 
hydrodynamic symmetries generates the infinite series of higher symmetries of the system 

(ED 

< = X(/- >v ),;>/;(./' + tvj(u) + Wj (u))). (5.61) 
i=i 

If R is a second order recursion operator, then all these symmetries are of the even order 
2N. 



6 Multi-component diagonal systems explicitly 
dependent on t or x 

6.1 Hydrodynamic symmetries of ^-dependent systems 

We consider the first order quasi-linear diagonal system of explicitly t-dependent equations 
with n > 3 

u\ = Vi{u,t)u l x , i = l,2,...,n (6.1) 

subject to the condition V{ ^ Vj for i ^ j. We search for hydrodynamic symmetries of this 
system with the Lie equations 

n 

< = £4M,zR z = l,2,...,n (6.2) 

i=i 

and x T = t T = 0. 

Define the functions 

Cij(u,t) = v iiU j(u,t)/(vj - Vi) (i^j), (6.3) 

rj,(tl)=Py(«,0). (6.4) 

Theorem 6.1 (see jl9], 20, 2^]). Diagonal n-component system (|6.1|) of the hydrody- 



namic-type with an explicit time dependence possesses an infinite set of hydrodynamic 
symmetries of the form (|6.2|) with a functional arbitrariness iff its coefficients meet the 
condition ( |5.12 ) and the condition 



[v iiU j(u,t)/(vi(u,t) -Vj(u,t))] t = /3v ijU j(u,t) (i ^ j) (6.5) 
with an arbitrary real constant (3. These symmetries are generated by the Lie equations 

vi r = A\u,t,x)vi x , i = l,2,...,n (6.6) 
with the coefficients Ai defined by the formulae 
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for 



and for (3 = 



t, x) = aj(u) exp < (3 



x + Vi(u, t)dt 



+ C 



Ai(u, t, x) = Oi(u) + C 



x + Vi(u, t)dt 



(6.7) 



(6.1 



where C is an arbitrary constant and the integrations with respect to t are performed at 
a constant value of u. Here the set of functions {cii(u)} is an arbitrary smooth solution of 
the linear system 

a i>uj (u) = r. j (w)(a i - a») (z 7^ j) (6.9) 

with the coefficients r*. defined by the formulae ( |6.3| ), ( |6.4|) . The solution manifold of 
the system ( |6.9|) depends upon n arbitrary functions Ci(u l ) of one variable which locally 
parametrize the set of hydrodynamic symmetries. 



6.2 Hydrodynamic symmetries of x-dependent systems 

Let coefficients of the diagonal system of (n > 3) equations explicitly depend on the 
coordinate x 

u\ — Vi(u, x)u x , i — l,2,...,n (6.10) 

and meet the condition V{ 7^ Vj (i 7^ j). 
Define the functions 



Cij(u, x) = v i)UJ (u, x)/ (vj - Vi) (i ^ j), 



(6.11) 

T^u) = cy(u, 0)vj(u, 0)/vi(u, 0). (6.12) 

Theorem 6.2 (see |T9],pD|,p2|). Diagonal n-component system ( |6.10| ) of the hydro- 
dynamic type with an explicit x-dependence possesses an infinite set of hydrodynamic 
symmetries of the form ( |6.2| ) with a functional arbitrariness iff its coefficients meet the 
condition (|5.12 ) with a substitution of Vi for Vi and the condition 



fa 1 (u,x)) uj /(v i 1 (u,x) -v^(u,x)) = Pfc 1 (u, x))^ (i^j) (6.13) 

with an arbitrary real constant (3. These symmetries are generated by the Lie equations 

u\ = Ai(u, t, x)u l x , z' = l,2, ...,n (6.14) 

with the coefficients A4 defined by the formulae 
for (3 ^ 



Ai(u, t, x) = Vi(u, x) < a,i(u) exp 



'> [ t + J v i (u, x)dx 



C 



(6.15) 
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and for (3 = 



Ai(u, t, x) = Vi(u, x) < di(u) + C 



t + 



u, x)dx 



(6.16) 



where C is an arbitrary constant and the integrations with respect to x are performed at 
a constant value of u. The set of functions {cii(u)} is an arbitrary smooth solution of the 



linear system (|6.9|) with the coefficients - defined by the formulae ( |6.11|) , ( |6.12|) . The 



extent of arbitrariness is the same as in the theorem 6.1. 



6.3 Invariant solutions and linearization 



Equations determining the invariant solutions of the systems ( |BTT1) and ( |6.10|) subject to 
the constraint u l x ^ for i = 1, 2, . . . , n are obtained from the Lie equations ( p.6| ) and 
( |6.14j ) with the invariance condition u\ = 0. They have the form Ai(u,t,x) = and 
Ai(u,t,x) = with the functions Aj and Ai defined by the formulae (|6~7|), ( |6.8| ) and 
( |6.15| ), (|6.16|) respectively where we put c = 1 without the loss of generality. Hence we 
obtain 

for the system (|6.1|) 



a,i(u) + exp < -f3 



x + 



(u, t)dt 



03 ^0), 



(6.17) 



ai(u) + x + / Vi(u, t)dt = {ft = 0), 



(6.18) 



and for the system ( |6.10| ) 



ai{u) + exp 



- '> [ t + J v i (u, x)dx 



(13^0), 



(6.19) 



aAu) + t + 



c 1 



(u,x)dx = (P = 0) 



(6.20) 



where i — 1, 2, . . . , n. 

Here the set of functions cii(u) is an arbitrary smooth solution of the linear system 
( |6.9D . Thus, the equations ( |6.17| )-( |6.20| ) determine a linearizing transformation for the 
systems (|6.1|), ( |6.10|) , i.e. the solution of these systems reduces to the solution of the 
linear system (|Q|). These equations determine the solutions u l = u l (x,t) of the original 
nonlinear system if the conditions of the implicit-function theorem are met. 

More complete results on the diagonal systems with explicit t or x dependence and 



the example of a new integrable system of this class can be found in the recent paper |22] 



51 



7 Conclusions 



The existence of an infinite-dimensional group of the hydrodynamic symmetries for the 
equations of the hydrodynamic type is an important property which provides the existence 
of linearizing transformations. The reason for this is that the degree of generality of the 
set of symmetries coincides with the degree of generality of the general solution set for 
these equations. Therefore, almost all solutions are the invariant solutions and they are 
obtained by standard formulae provided the symmetries are already determined. Such a 
formula gives a linearizing transformation reducing the original non-linear problem to the 
linear problem of determining the symmetries. The additional property is the existence 
of the recursion operator which makes it possible to solve partially the linear problem by 
constructing infinite series of its solutions and hence solutions of the original non-linear 
equations. The existence of such an operator also has a group-theoretical basis, since 
the recursion operator is completely determined by the symmetries of the determining 
equations for the hydrodynamic symmetries, i.e. by the 'symmetries of symmetries". 

This shows a group-theoretical origin of linearizing transformations and of the integra- 
bility property by which we mean a possibility to construct infinitely many exact solutions. 
The Hamiltonian structure, if it exists, does not improve the integrability properties of 
the equations of the hydrodynamic type. We have to conclude that the symmetry is the 
major necessary property that insures the integrability which was the original idea of 
S. Lie. 
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